Copy  No.      80 


NEW  YORK  UNIVERSITY 

IN  ^  1  OF  MATHEMATICAL  SCIENCES 

LIBRARY 

25  Waverly  Place,  New  Yofk  3,  N.  Y. 


AMP  Report  137.1R 
AMG-NYUNo.  128 


THEORETICAL  STUDIES  CONCERNING  THE  HVDROPULSE 

(IDEAL  MECHANICAL  PERFORMANCE  CHARACTERISTICS) 


WITH  THE  APPROVAL  OF  THE  OFFICE  OF  THE 

CHAIRMAN  OF  THE  NATIONAL  Dr.FcN' E  RES'iARCH 
COMMITTdE,  JHiS  REPORT  IIAS  SEEM  DECLAcSIFIED 
BY  THE  OFFICE  OF  SCIENTIFIC  RESEARCH  AND 
DEVELOPMENT. 


Prepared  for  the 

APPLIED  MATHEMATICS  PANEL 

NATIONAL  DEFENSE  RESEARCH  COMMITTEE 

By  the 

Applied  Mathematics  Group 

New  York  University 


X 


This  document  contains  Information 
affecting  the  National  Defense  of  the 
United  States  within  the  meaning  of  the 
Espionage  Act,  U.  S.  C.  50;  3!  and  32. 
Its  transmission  or  the  revelation  of  its 
contents  In  any  manner  to  an  unauthorized 
person   is   prohibited    by  law. 


■eo-N-PRffliPPfrA  t 


aJ^^ 


^' 


July  1945 


Prepared  in  oormeotion 
with  Project  NA-195 
under  Contract  OEMsr-946 


:AL 


ApproTed  for  Distribution 
l&rren  Weaver 
Chief,  Applied  Uathematios  Panel 


Distribution  List 
AMP  Report  137.1R 


Copy  No< 
1-7 
8-34 
35 

S6 

57  -  40 


41  -  45 


46  -  47 


48-54 


55  -  59 


60  -  61 


62 


Office  of  Executive  Secretary,  OSRD 

Liaison  Office,  OSRD 

Office  of  Research  and  Inventions 
Att:   Lt.  J.  H.  Wakelln 

Commanding  General,  AAP 
Attt    T,  von  Karman 

Chief,  Bureau  of  Asronautica 
1  Comdr.   J.  S.  Warfel 
1  Lt,  F.  Ao  Parker 
1  B,  S,   Roberts 

Bureau  of  Aeronautics  Representative  in  Los  Angeles 
Att:  Aerojet  Coirp, 

U.S.  Naval  Engineering  Experimental  Station,  Annapolis 
1  Lt,  Patton 
1  Ensign  L,  B,  EdeLnan 

Chief,  Bureau  of  Ordnance 
1  R,   S,  Burington 
1  R,  J.  Seeger 

1  Comdr.  Levering  Smith,  Re2d 
1  Comdr.  W,  A,  Walter,   Re3 
1  Capt,  D,  P,  Tucker  Re4g 
1  Capt,  E.  U.  Crouch  R06 

Conraanding  General,  AAF  Materiel  Coinnand,  Weight  Field 
1  Lt,   Col.  P.  F.  Hay 
1  Major  J.  P.  AuWerter 
1  Capt,  J.  Healy,  Armament  Laboratory 
1  Lt,  R,  C,  Bogert 
1  Bell  Aircraft  Corp, 

J,  E,  Jaokson,  Chairman's  Office,  NDRC 
1  Busblebee  Project 

Guided  Missiles  Committee,  JCS 


00. 


"UL 


CONFIDENTIAL 


63  -  66    Colonel  L.  A,  Skinner,  Ordnance  Department  Representative,  C.I.T, 
1  C.  Millikan 
1  L,  G,  Dunn 
1  J.  V.  Charyk 

67  -  69    Office  of  the  Chief  of  Ordnance 
1  Colonel  G.  W.  Trichel 
1  Major  J,  F,  Miller,  Jr. 
1  H.  M.  Morse 

70  Commanding  General,  AAF  Proring  Ground  Conanand,  Eglin  Field,  Fla, 

Att:  Brig,  Gen,  Grandison  Gardner 

71  -  72    Ordnance  Research  Center,  Aberdeen  Proving  Ground 

1  Technical  Library 
1  R.  J.  Walker 

73  -  74    F.  L,  Hovde,  Chief,  Division  3 
IE.  C,  Watson 
1  J,  B,  Rosser 

75  J,  T,  Tate,  Chief,  Division  6 

76  W.  R,  Kimer,  Chief,  Division  9 

Atti  J,  W,  Williams 

77  W,  Weaver,  Chief,  Applied  Mathematics  Panel 

78  T,  C,  Fry,  Acting  Chief,  Applied  Mathematics  Panel 

79  -  80    R,  Courant  '^ 

81  0.  Veblen 

82  Q.  C.  Evans 
85  S.  3.  Wilks 
84  M.  Rees 


CONFIDENTIAL 


CONFIDENTIAL 


PREFACE 


The  present  report,  by  the  New  York  University 
Group  of  the  Applied  Mathematics  Panel,  was  tindertaken 
in  connection  with  Project  No.  NA-195  (Bxireau  of 
Aeronautics,  Navy  Department).   It  was  carried  out, 
under  the  general  guidance  and  editorship  of  Dr. J. J. Stoker, 
by  M.  Shifftnan,  D.  C.  Spencer,  B.  Friedman,  E.  Bromberg 
and  E.  Isaacson  in  consultation  with  Dr.  P.  Zwlcky  of  the 
Aerojet  Corporation. 

The  report  contains  a  mathematical  analysis  of  a 
variety  of  problems  concerning  the  mechanics  of  the  exhaust 
stroke  in  a  hydropulse  motor.  A  broad  and  basic  approach 
to  the  problems  was  motivated  to  a  considerable  extent  by 
the  desire  to  avoid  too  specific  assxmiptlons  with  regard 
to  ^he  engineering  aspects  of  the  hydropulse  development. 


R.  Courant 

Director  of  Research 

Contract  0EM3r-945 
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Ideal  Mochanlcal  Performance   Gharaeterlatlea 

of  the   qydropulae 
Section  lo     Introduction  and   aummary. 

In  this  report  ^e  are   interested  in  certain  aspects 
of   the  problems  which  arise   in  dealing  with  jet  propulsion 
devices   to  be  used  in  water  —  in  other  words,  with  the 
counterparts   in  water  of  the    jet  devices   in  air,   such  as 
rockets,    the   athodyd,    the  buzz  bomb  motor,   etc.      The 
manner  in  which  such  devices  function  can  be  described  in 
a  general  way  as  follows:      The  propulsive  element  is   a  duct 


V 


TUBE 


Valves 


Figure   1.1 

open  at  both  ends   (as   indicated  la  Figure   1.1).     The  duct 
is  assumed  to  be   traveling  through  the  water  with  a   certain 
velocity  V  (presumably  it  is  propelling  a  ship,   a  torpedo, 
or  other  craft).     At  a   certain  Instant  a   chemical  is 
Injected  into   the  duct  which  reacts  with  the  water  to  produce 
a  volume  of  gas,   as  indicated  In   the  flgvire.      The  gas,  on 
expanding.   Imparts  momentum  to  the  water  and  a  thrust  in  tha 
forward  direction  will  result  if,   for  example,    the  part  of 
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the  duct  in  front  of   the  gas   la   sealed   off  by  closing  valves. 
If  the  forward  velocity  is  high  enough  and  if  the   intake 
part  of  the  duct  is  a  properly  designed  expanding  nozzle, 
valves  may  not  be  needed  to  insure  a  forward  thrust  and  the 
fuel  could  be   injected   continuously.      Such  a  device  would 
correspond  to  the   athodyd  in  air.     However,  at  moderate 
velocities  V,    it  is   likely  that  valves  would  be   necessary 
for  efficient  operation  of  such  a  duct  device.      Of  course, 
sustained  operation  of  the  device  with  valves  requires 
the  valves   to  open  at  certain   times  so  that  the  expansion 
chamber  can  be  refilled  with  water  after  the  gas  has  expanded 
and  pushed  a   predetermined  amount  of  water  out   of  the   tube. 
Once   the    "expansion  chamber"    has  been  refilled,    the  valves 
would  be   closed,   fuel  would  be   injected,  and  the  cycle 
would  repeat. 

In  what  follows  we  assume  that  valves  will  be  provided, 
but  we  do  not  assxime    that  we  know  when  and  how   they  function  - 
whether  they  open  automatically  when  the  gas  pressure  reaches 
a  certain  value    (the  hydro-resonator),   or  are  controlled  by 
some   Independent  mechanism  ( hydropulae ) .        Our  investigations 
have   therefore  been  confined  for  the  most  part   to   those 
phenomena  which  can  be  treated  without  knowing  in  detail 


This   terminology  we  have   taken  from  the  report  of  P. 
Zwicky:     Remarks   on   the  Basic  Theory  of   the   Qydropulse   as 
a  Propulsive  Power  Generator,     Aerojet  Report  No.  R-41 
(Conf.).     In  fact,  most  of   the  work  to  be  discussed  In 
what  follows  was  undertaken  at  the   suggestion  of  P.   Zwicky 
and  in   consultation  with  him.     Two  other  reports  of  Interest 
for  the  mechanical  performance  of  the  hydropulae  have 
appeared:      1)   The   Solution  of  the  Differential  Equation 
for  the  hydropulae    (Conf.),   Aerojet  Report  No.   R«-43,  by 
K.  Poelsch  and  Knox  T.  Mlllsaps,   2)   A  Preliminary  Study 
of  the  5ydropulse    (Conf.),   Jet  Propulsion  Note  No.   19, 
Bureau  of  Aeronautics  Project,  U.S.   Naval  Eng.  Exp. 
Station,  Annapolis,  Md.,  written  by  Ensign  L.   B.  Sdelman. 
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just  how  the  complete  cycle  of  operations  will  be  carried  out. 
Our  calculations  should  be  helpful  in  making  decisions  on 
this  point.  Roughly  speaking,  what  we  have  done  is  the 
following:  We  have  ignored  the  part  of  the  cycle  concerned 
with  the  refilling  of  the  "expansion  chamber",  (except  in 
Section  12)  and  concentrated  our  attention  on  that  part  of 
the  cycle  during  which  the  gas  is  expanding  and  pushing  water 
out  of  the  tube.  We  consider  the  motion  up  to  tie  time  when 
the  gas  has  expanded  down  to  the  exit  pressure,  and  define 
as  the  time  of  stroke  —  or  simply  the  stroke  —  the  time 
required  for  such  an  expansion.  In  carrying  out  the  calcula- 
tions certain  simplifying  assumptions  are  made  which  are 
discussed  in  Section  2,  immediately  following  this  intro- 
duction. Briefly  speaking,  these  assumptions  lead  to  a 
"hydra laic  theory",  in  which  the  water  is  treated  as 
incompressible  and  the  velocity  distribution  over  any  cross- 
section  of  the  tube  is  replaced  by  the  average  over  the 
section.  In  addition,  the  gas  is  assumed  to  be  uniform  in 
state  throughout  its  volume  at  each  instant  of  time,  and  to 
expand  adiabatically  in  such  a  way  that  the  interface  between 
gas  and  water  is  at  all  times  a  plane  surface  at  right  angles 
to  the  axis  of  the  tube.   The  pressure  at  the  exit  of  the 
tube  is  assumed  to  be  constant  throughout  the  stroke.   On  the 
basis  of  these  assumptions  we  study  the  motion  of  the  water 
in  the  tube,  the  impulse  delivered  over  the  stroke,  and  the 
time  of  stroke  —  that  is,  we  study  the  mechanics  of  the 
hydropulse  operation,  but  leave  aside  the  chemical  and  thermo- 
dynamic problems  concerned  with  the  production  of  the  gas* 
In  order  to  calculate  efficiencies  we  simply  assume  that  a 
certain  mass  of  a  perfect  gas  is  present.   These  assumptions 
which  seem  rather  reasonable  for  the  purpose  of  a  preliminary 
attack,  should  yield  the  ideal  mechanical  performance  data 
for  devices  of  this  type. 
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Wo  outline   the   contents  of   this  report  section  by 
section  together  with  the  main  conclusions  reached  In  eachs 

(Section  2*      Contents  discussed  above.) 

Section  5«     Formulation  of  the  hydraulic   theory*      The  equation 
of  motion  of  the  water  during  a   single   stroke  Is  derived  from 
the  relation  for  the  energy  balance.      The  formula  for  the 
Impulse   Is  also  obtained.      (This   is  essentially  the   same  as 
the  theoiT"  derived  in  the  previously  quoted  report  of  P. 
Zwicky. )     The   theory  leads  to  an  initial  value  problem  for 
a  non-linear  ordinary  differential  equation  of  second  order. 

Section  4.     Introduction  of  dimensionlesa  variables  and 
parameters.      The  essential  variables  and  parameters  are 
expressed  in  dimenslonless  form.     In  all  cases  Involving  thQ 
determination  of  quantities  as  functions  of  the  time   through 
actual  integration  of  the  equations  of  motion   (as   in  part  of 
Section  7  and  in  Sections  9  to   12  incl.)   these  variables  ai^ 
very  convenient. 

Section  5.     Duct  of  variable  cross-section i     Approximate 
formulas.      Shape  of  duct  for  maximum  impulse.     Although  the 
hydraulic   theory  as  developed  in  Section  3  is  quite  elementary 
from  a  mathematical  point  of  view,   the  amount  of  calculation 
necessary  for  a  complete  discussion  of  the  motion  of  the 
water,   the  impulse  I     per  stroke,   the  ratio  of  impulse  jjer 
stroke   to  mass  of  fuel  injected  per  stroke    (the  "effective 
velocity"  or  "specific  Impulse"   Ig/^g)#    the  time  of  stroke   t   , 
and  the  average  thrust  F     =  I_/t_  in  their  dependence  on  all 
of  the  essential  parameters  would  be   truly  enormous.      The 
reason  for  this  is   that  the  equation  of  motion  can  not  be 
Integrated  explicitly,   so  that  numerical  Integration  methods 
must  be  used.      However,   it  turns  out  that  the  quantities   just 
now  enumerated  above,  which  are   the  most  important  ones  for 
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practical  purposes,  can  be  expressed  by  approximate 

formulas  which  do  not  involve  a  complete  integration  of  the 

equation  of  motion.  Section  5  contains  a  derivation  of  these 

formulas,  together  with  proofs  that  the  formulas  yield 

rigorous  upper  bounds  for  I_  and  f  and  a  rigorous  lower 

s  s 

bound  for  t„. 
s 

By  comparison  with  the  results  obtained   through 

accurate  numerical  Integration  of  the  equation  of. motion 

in  a  large  number  of  cases   the  formulas  are  known  to  furnish 

results  which  are  accurate  within  2  per  cent  for  I     and 

a 
5-15  per  cent  for  t^  as   long  as   the   initial  pressure  ratio 

Is  greater  than   10  and  for  ducts  which  do  not  depart  in 
shape  too  much  from  that  of  the   straight  tube  of  uniform 
cross-section. 

Section   5  also   contains   a   discussion  of  the   question 
of  maximizing   the   Impulse  with  respect  to   changes  in  shape 
of   the  hydropulse   tube.      In  case   the  forward  velocity  of 
the  hydropulse   is   small,   it  Is  not  difficult   to  show   that 
the   straight   cylindrical  tube  Is  as  good  as  any  other  for 
this  purpose.      This  conclusion  follows  from  the  fact   that 
our  approximate  formula  for  I     yields  an  upper  bound  for  all 
shapes  of  tubes  and   the  fact   (determined  by  actual  calcula- 
tion) that  the  upper  bound  is  attained  within  a  per  cent  or 
two  In  the   case  of  the    straight  tube.      Bowever,   If  the 
forward  speed  of   the   hydropulse   is  high,   so   that  the  velocity 
of   the  water  in  the   tube   Is  high  at  the  beginning  of   the 
stroke,   the   Impulse  per  stroke  I^   Is  affected  somewhat  by 
changes  In  the  shape  of  the   tube.     We  have  not  obtained  a 
final  result   In  these   cases,  but  some   conclusions  from  the 
present  section,   combined  with  a  few  numerical  results 
from  Section  8  on   tubes  with  conical  sections.   Indicate   that 
the  best  shape  of  tube  for  maximum  Impulse  during  the  power 
stroke  at  a  given  forward  speed  would  be  one  with  a  flared-out 
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forward   part   joined   to  a   straight   cylinder  extending   to   the 
exit,   as   Indicated    In   the   accompanying   sketch. 

Dt  recti  on    or     motion 


Vol 


vea 


Section  6«      The   tube  of  constant  cross -section:     Formulas 
showing  dependence  of  essential  quantities  on  all  parameters* 
The  main  bulk  of  our  calculations  concern  the   tube  of  constant 
cross-section.      Our  objective   is  a   complete  discussion  of  the 
behavior  of  all  Important  quantities  when  any  or  all  parameters 
are  varied.     We  have  already  stated   that  approximate  formulas 
are  available   for  tubes  of  any  shape.      These  formulas  are 
developed  explicitly  in  Section  6  for  the   special  case  of 
the   straight  tube   In  such  a  way  as   to  make  a  complete  dis- 
cussion possible  for  the  full  ranges   of  all  essential 
parameters.     In  particular,  formulas  for  the  Impulse  I   , 
the  effective  velocity  Ig/*^-*    (which  is  a  measure  of  the 
efficiency  of  the  device),   the   time  of  stroke  t„,  and  the 
average   thrust  P     =  I_/t_  are  obtained.      These  formulas 
exhibit  explicitly  the  dependence  on   the  following  parameters: 
1)   the  Initial  pressure  ratio  P^/Pq*  2)   the  ratio  M  yM     of 
the  mass  of  gas  M     Injected  per  stroke   to   the  initial  mass 
JL  of  water  in  the   tube,  3)   the   initial  velocity  U.   of  the 
water  in  the   tube,     4)   the  adlabatle  constant  7,5)   the 

The  initial  velocity  TJ^  of  the  water  in  the  tube  would  be 
approximately  equal  to  the  forward  velocity  of  the  hydro- 
pulse. 
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Initial  temperature   and  molecular  weight  of   the  gas. 

At  the  end  of  Section  6  graphs  and  a  table  are  given 
which  compare   the  results  obtained  by  the  approximate 
formulas  with  those  obtained  by  accurate  numerical  integration 
of  the  equation  of  motion.      The  approximate  formulas  are 
found  to  be  amply  accurate  for  our  purposes. 

Section  7>      The   tube  of  constant  cross-section.      Performance 
curves,  with  discussion.     This   section  is   perhaps   the  most 
important  one   in  the  report.      The  behavior  of  the  essential 
quantities  with  respect  to  variations  of   the  parameters  is 
given  in  the  form  of  graphs  placed  at  the  end  of  the   section. 
A  brief  summary  of   the  main  results  follows: 

a)     The  effective  velocity  Is^g,   or  ratio  of   impulse 
per  stroke   to  mass   of  fuel  injected  per  stroke.      The  effective 
velocity  increases  quite  slowly  once   the   initial  pressure 
ratio  Pi/Pq  becomes  greater  than  about  30,   particularly  for 
the  higher  initial  velocities.      Thus   there  would  be  no  great 
advantage  from  the  point  of  view  of  fuel  conservation  in 
going  to  pressure  ratios  much  above  50.      The  effective 
velocity  decreases   steadily   to  zero  as  M     increases   steadily 
to  infinity.      On  the   other  hand,  as  M     tends  to   zero,   the 
effective  velocity  increases   steadily  and  tends   to  a  definite 
limit  which  is  finite   if  the   initial  velocity  JJ^  is  not  zero, 
but  infinite   if  U.    is   zero.      It  is  therefore  clear  that 
efficient  operation  is  promoted  by  keeping  the  amount  of  fuel 
injected  per  stroke  as   low  as   that  which  is  compatible  with 
the  attainment  of  other  objectives.     It  is   clear  that  M     must 
be   large  enough  so   that  a  sufficiently  high  thrust  can  be 
attained.     As  the   initial  velocity  TJ^  increases,   Ig/^g  decreases 
and  has   the   limit  zero.     This  is  to  be   expected,   since   the 
thrust  augmentation  is   lowered  with  increase  in   the   initial 
velocity  of  the  water. 

The  actual  numerical  values  for  the  effective  velocity 
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in  practically  reasonable   cases  are  of  the  order  of  several 
hundred  thousand  feet  per  second.      These  are  high  values 
compared  with  those  achieved  with  jet  devices  In  air,   the 
difference  being  due  to   the  very  great  density  of  water  as 
compared  with  air. 

b)  The   time  of  stroke  tsJ      The  time  of  stroke   is  not 
very  sensitive   to  changes  in  the   initial  pressure  ratio  P^/Pe 

if   this   ratio  exceeds   50  or  40«      When   the  fuel  ratio  M  /M 

.^.»— — — — — — — — — — — — — — —  g/  w 

alone  is  varied  the  following  behavior  of  t_  is  found j  As 

s 

M  M     tends  to  zero  so  also  does  the  time  of  stroke  tg,  while 

°  /M 

t     tends   to   infinity  like  /«S.  when  M  /M     tends   to  infinity. 

S  V   'vj         o 

That  t  should  tend  to  zero  with  M  seems  at  first  sight 
s  g 

rather  strange,  but  it  is  readily  explained  by  the  fact  that 

the  initial  pressure  ratio  is  assumed  to  be  held  fixed  while 

M  tends  to  zero.  This  would  be  difficult  to  achieve  in 
g 
actual  practice. 

c)  The  average  thrust  P^  =  13/tg.  The  average  thrust 

increases  steadily  with  increase  in  the  initial  pressure 

ratio,  as  we  would  expect.  When  the  fuel  ratio  M  /M^  alone 

is  varied,  the  average  thrust  P_  behaves  as  follows:  P 

M 

decreases  steadily  from  a  finite  value  for  jp*  -  °°  and  tends 

w 

as  M  /M  -»•  0  to  another  finite  value  which  is  half  the  value 
g/  w 

for  M  M.     -  00  ,  Again  this  somewhat  paradoxical  behavior  of 

F.  as  the  quantity  of  fuel  M  injected  per  stroke  approaches 
s  g 

zero  is  explained  by  the  fact  that  the  Initial  pressure  ratio 

is  assumed  to  be  maintained  as  M^  -^  0,  and  this  could  not  be 

g 
achieved  in  practice.  Nevertheless,  in  view  of  the  fact  that 

the  efficiency  as  measured  by  the  effective  velocity  I./^I- 

has  its  maximum  for  M  yM  =  0,  the  following  important 

conclusion  is  strongly  indicated:   The  hydropulse  should  be 

operated  at  high  frequency  with  injection  of  small  amounts 
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of  fuel   per  stroke.      The   quantity  of  fuel  In.jected  pr  stroke 
should  be  kept  as   lov?  as  Is   compatible  with  maintenance  of  a 
minimum  Initial  pressure  ratio   (of  the  order  of  30  or  40  If 
possible)   30  that   the  average   thrust  will  be  maintained  at 
the  value  necessary  to  overcome   the  drag  on   the  hydropulse* 

d)  Gas  pressure  and  exit  velocity  of  the  water  aa 
functions  of  the   time.     In  this  part  of  Section  7  we  report 

a  few  of  our  calculations  on   the  behavior  of  the  gas  pressiire 
and  velocity  of  the  water  as  functions  of  the   tlmfe.      These 
results  are  of  Interest  In   themselves.      They  also   lead  to   the 
development  of  the   approximate  formulas.      The   principal  ob- 
servation to  be  made   is   that  the  gas   pressure  decreases  very 
rapidly  Initially,  while   the  velocity  of   the  water  Increases 
very  rapidly  at  the  beginning.     Both  quantities   then  remain 
nearly  constant  after  the  first  short  interval  of   time  at  the 
beginning  of  the  stroke.     A  number  of  graphs  exhibiting   this 
behavior  in  a  few   typical  cases  are  given  at  the  end  of 
Section  7.     In  Sections  9  and  10  two  possible  means  of  smooth- 
ing out  these   initial  abrupt  variations  in  pressure  and  in 
velocity  of  the  water  are  discussed. 

e)  Effect  of  linear  buildup  (in  time)  of  the  gaa 
pressure.     In  all  but  this   portion  of  Section  7  we  assume 
that  the   initial  high  pressure  of  the  gas  la   created  instan- 
taneously.    Actually  the  gas   pressure  would  probably  be 
built  up  gradually  to  its  maximum  after  injection  of  the  fuel. 
In  what  manner  the  pressure  would  rise   in  actual  practice 

is  not  known  —  it  would  depend  on  the  chemical  reactions  and 
the  thermodynamic a  of  the  process.  In  this  part  of  Section  7 
we  give  the  results  of  a  few  calculations  assuming  the  gaa 
pressure  to  rise  linearly  in  time  to  its  maximiun,  after  which 
the  gaa  expand a  adiabatically  down  to  exit  pre a  sure  (as  in 
all  of  our  other  calculationa ).  For  the  caae  of  zero  initial 
velocity  of  the  water  in   the   tube   curvea  are  given  from  which 
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the  gas  volume  and  rate  of  expansion  are  determined  for  any 
given  rate  of  increase  in  pressure.  The  time  required  for 
the  pressure  to  reach  its  maximum  was  taken  as  about  30  per 
cent  of  the  entire  stroke  in  one  numerical  case  and  40  per 
cent  in  another.  The  time  of  stroke  and  impulse  per  stroke 
were  computed  in  a  number  of  cases  and  the  results  compared 
with  cases  in  which  the  gas  pressure  reaches  its  maximtim 
instantaneously.   The  results,  though  inconclusive,  are 
interesting.   In  these  few  cases  we  find  a  gradual  buildup 
in  pressure  to  be  favorable  for  the  operation  of  the  hydro - 
pulse,  since  the  impulse  delivered  per  stroke  for  a  given 
quantity  of  fuel  appears  to  be  substantially  larger  than  in 
the  case  of  instantaneous  rise  of  pressure.   (However,  our 
method  of  comparing  the  quantities  of  fuel  may  be  open  to 
serious  objection  on  physical  grounds.  For  details,  see 
Section  7e). 

If  the  last  conclusion  were  to  be  found  generally 
valid,  it  would  mean  that  our  efficiencies  (which  were  all 
computed  on  the  assumption  of  instantaneous  pressure  rise) 
are  on  the  conservative  side. 

Section  8.   Tubes  with  conical  portions.   The  results  of 
calculations  giving  the  Impulse  and  time  of  stroke  for  three 
different  types  of  tubes  having  conical  portions  In  front 
are  presented  and  compared  with  those  for  the  straight 
cylindrical  tube.  The  purpose  and  the  main  results  of  this 
section  have  already  been  stated  in  the  course  of  our  remarks 
above  referring  to  Section  5. 

Section  9.   Straight  tube  with  mass  and  spring  for  elastic 
energy  storage.  One  of  the  results  of  the  calculations  of 
Section  7  is  that  the  gas  pressure  (and  consequently  the 
thrust)  and  the  velocity  of  the  water  in  the  tube  change  very 
rapidly  at  the  beginning  of  the  exhaust  stroke.  In  order  to 
promote  efficiency  it  seems  likely  that  a  more  even  distribu- 
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tion  of  th9  thrust  and  of  the  exit  velocity  of  the  water 
would  be  desirable.   One  possible  way  to  do  this  would  be  to 
provide  for  an  elastic  storage  of  energy  during  the  early 
part  of  the  stroke  which  would  then  be  delivered  again  in  the 
later  part  of  the  stroke.  Just  what  means  should  be  taken 
to  accomplish  this  in  practice  is  far  from  clear.  What  we 
have  done  is  to  investigate  a  mechanical  model  in  order  to 
find  out  if  possible  what  can  be  accomplished  by  such  means. 
The  model  Is  obtained  by  assuming  that  the  left  end  of  the 
exhaust  chamber  is  closed  off  by  a  piston  of  a  certain  mass 
which  is  attached  to  a  spring  (instead  of  a  rigid  wall  at  the 
left  end.)  The  calculations  were  carried  out  for  various 
values  of  the  mass.  It  turns  out  that  the  time  distribution 
of  the  thrust  and  the  velocity  of  the  water  can  be  very 
considerably  influenced  through  proper  choice  of  a  mass  and 
spring.  In  particular,  the  abrupt  initial  variations  in  these 
quantities  can  be  smoothed  out  rather  successfully.   (It 
might  be  noted  that  a  spring  acting  alone  without  a  mass  is 
not  effective. ) 

It  is  found  that  the  operation  of  the  hydropulse  might 
be  rather  erratic  for  certain  values  of  the  parameters  since 
the  time  of  stroke  may  be  discontinuous  for  special  values 
of  the  parameters. 

Section  10.   Spring  effect  obtained  from  a  layer  of  gas.  One 
method  of  approach  to  the  problem  of  providing  elastic 
storage  of  energy  during  part  of  the  stroke  might  be  to  inject 
the  fuel  in  such,  a  way  as  to  create  two  gas  volumes  separated 
by  a  mass  of  water.  Various  cases  of  this  kind  are  treated 
in  Section  10.  The  number  of  parameters  which  could  be  varied 
in  these  eases  is  rather  large,  and  our  calculations  cover 
only  a  small  number  of  the  possibilities.   However,  they  do 
show  that  such  a  method  could  be  used  to  modify  to  some  extent 
the  time  distribution  of  thrust  and  velocity. 
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Section   11 »      The  valve lea a   hydropulae.     Except  In   this   and 

the   section   to  follow  we   assume  always   that   the   left  end  of 

the   hydropulse    tube   is   closed  by  valves  of   some  kind.      In 

Section   11  the   possibility  of  operating  a   hydropulse   without 

valves   is   considered.      The   tube   is  assumed   to  be   open  in 

front  but  to  have  a  constricted  front  portion.     At  the   instant 

of  injection  of   the  gas   at  relatively  high  pressure    (behind 

the   constriction,   of  course)  water  is   assumed   to  be  flowing 

through  tl^   tube.      The  motion  of  the  water  in  the  tube  is 

analyzed,  on  the  basis  of  our  usual  assumptions,  up  to  the 

time  when  the   gas   pressure  vhas  decreased   to    the  value   at  the 

exit.      The   solution  of   the   problem  for  any  one    set  of 

'it- 
parameter  values  requires   the  numerical  integration     of  a 

pair  of  second  order  non-linear  differential  equations. 

Pour  cases  were   treated.      The  general  conclusion  reached  is 

the  following  J     As  far  as   the  total  impulse  for  a  single   stroke 

is   concerned,   a  hydropulse  without  valves   could  be  designed 

so  that  its  performance  would  compare  favorably  with  that  of 

the  hydropulse  with  valves.      To  do   this,  however,  a  fairly 

long  narrow  forward  or  nose   section  would  seem  to  be  needed. 

(In  our  calculations  we   assumed   the   cross-section  area  of  the 

nose   to  be  one -fourth  that  of  the  main  part  of  the  hydropulse 

tube.)     As  we   shall  see   shortly,   such  a  constricted  forward 

portion  of  the   tube  would  be  very  unfavorable  for  the 

scavenging  or  refill  portion  of  the  whole   cycle  of  operations. 

Section  12.     Scavenging.     We  quote  from  Section   12 J      "In  this 
report  we  have  in  general  concentrated  our  attention  on  a 
single  exhaust  stroke  of  the  hydropulse  and  have   left  aside 
problems   concerning   other  phases   of   the   full  cycle,    in  particu- 
lar,  the  phase  of    'scavenging'.     By  this  la  meant  the 
portion  of  the   cycle  during  which  the   tube   la  being  refilled 


The  method  used   is  explained   in  Appendix  I. 
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with  water  in  preparation  for  the  injection  of  fuel  and  the 
beginning  of  a  new  exhaust  stroke."  —  "The  problem  we  wish 
to  consider  is  the  following:  We  assume  that  the  exhaust 
stroke  has  just  been  completed,  i.e.  that  the  gas  has 
expanded  (with  the  valves  closed)  until  its  pressure  has 
decreased  to  the  value  p.  on  the  outside  of  the  tube  and  the 
water  in  the  exhaust  chamber  has  acquired  a  certain  velocity. 
An  instant  later  the  gas  will  be  at  a  pressure  less  than  p 
(since  the  inertia  of  the  moving  mass  of  water  will  lead  to 
a  further  expansion  of  the  gas,)  so  that  valves  could  be 
opened  automatically  by  the  pressure  of  the  water  in  the 
forward  part  of  the  tube.'  Our  purpose  in  this  section  is  to 
study  the  subsequent  motion  of  both  layers  of  water,  on  the 
basis  of  our  usual  assumptions  regarding  the  behavior  of  gas 
and  water  layers." 

Again  the  solution  of  the  problem  requires  the  numerical 
integration  of  a  pair  of  second  order  differential  equations 
involving  a  considerable  number  of  parameters.  We  have 
carried  out  the  calculations  (by  the  methods  of  Appendix  I) 
in  one  case,  using  what  seem  to  be  reasonable  values  for  the 
parameters.  It  turns  out  that  the  scavenging  operation  la, 
generally  speaking,  very  slow.  We  find  in  our  nvnnerical  case 
that  the  time  required  for  the  gas  pressure  to  reach  its 
minimum  after  opening  the  valves  is  three  times  as  long  as  the 
time  for  the  previous  power  stroke,  and  even  then  the  refilled 
volume  would  be  less  than  that  of  the  gas  layer  at  the  instant 
of  opening  the  valves. 

For  efficient  operation  of  the  hydropulse,  it  is  clear 
that  the  scavenging  phase  should  not  take  up  too  large  a 

This  Implies,  of  course,  that  the  tube  extends  a  certain 
distance  in  front  of  the  valves.  The  tube  is  also  taken 
to  be  a  cylinder  of  uniform  section  throughout. 
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fraction  of   the   full  cycle.      Our   calculations,   although  they 
cover  only  one   case,    thus    indicate    that  the   scavenging   phase 
In  the  operation  of  the   hydropulae  may  present  problems  of 
crucial  importance.      The   latter  part  of  Section  12  discusses 
a  few  possibilities  for  overcoming   the  difficulties,   which  in 
their  turn  lead   to  new  problems    (more   careful  consideration 
of   the   actual  flow   into   the  mouth  of   the   tube,   for  e  :x  ample ) 
which  will  require   investigation. 

Section  15.      Effect  of  finely  distributed   gas  bubbles.      In 
practice,   the  fuel  is   likely  to  be   injected  into  the  water  in 
such  a  way  as   to  result  in  the   creation  of  a  large  number  of 
minute  bubbles  distributed  over  a   certain  volume  of   the  water, 
rather  than  in  one  single   layer,   as  was  assumed   in  earlier 
sections.     The  motion  of  the  water  can  also  be   treated  without 
difficulty  in  this  case.      This  problem  is   treated  in  Section 
13  on  the  basis  of  a  number  of  simplifying  assumptions.     The 
main   conclusion  reached   is   that   the  motion   of  the  water  in 
the   tube  is   the  same  in   this   case  as   it  would  be  for   the 
corresponding  case  of  a  single   layer  of  gas   in  a   tube  of 
smaller  total  volume.      In  fact,    the  volume   of   the    tube  with 
one   layer  of  gas   should  be  diminished  by  one -third  of   the 
volume  of  that  part  of   the  water  ishich  contains   the  gas 
bubbles  in  order  to  make   the  two   cases   identical  with  regard 
to   total  impulse  delivered  per  stroke,   etc.      In  other  words, 
only  two-thirds  of   the   volume  of    the  water   containing   the  gas 
bubbles  is  effective  for  thrust  augmentation. 

Appendix  I.  The  results  of  this  report  entailed  a  great  deal 
of  nviraerical  computation  of  various  kinds  about  which  nothing 
much  is  said  in  the  main  part  of  the  report  —  to  have  done 
so  in  any  detail  would  have  at  least  doubled  the  length  of  an 
already  overlong  report.  However,  a  few  remarks  on  numerical 
methods  of  two  different  kinds  have  been  placed  in  two 
appendices.      Appendix  I  explains   the  method  of  finite  differ- 
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ences  which  was  used  for  the  ntmerlcal  integrations  of  the 
systems  of  differential  equations  encountered  In  Sections  9 
to  12  Inclusive*      This  method  is  not  at  all  a  ne«  or  unknown 
one,  but  it  was  felt  that  the  future  usefulness  of   this 
report  might  be  enhanced  by  the   inclusion   of  a  method  which 
we  found  convenient  and  efficient  for  solving  this   type  of 
problem.     Putxire  developments  might  well  require  refinements 
in  the   treatment  of  the   hydro pulse  which  could   lead   to  still 
more   complicated  aya terns  of  differential  equations   than  those 
in  Sections   9  to   12.      However,   the  method  explained  in 
Appendix  I  la   capable   of  dealing  with  such  cases  by  general- 
ization in  a  rather  obvious  manner. 

Appendix  II «     For  dealing  with  simpler  cases   in  which  only 
one  differential  equation  occurs   the  method  of  finite 
differences   can  be  replaced  by  a  method  which  gives   accurate 
results  with  very  much  less  labor  in  computation.      This 
method  is  explained  in  Appendix  II  for  the  special  ease  of 
the   straight  tube.     It  is,  however,  equally  useful  for  tubes 
of  varying   cross-section,   and   this  may  be   of  some  use   in   the 
future. 


CONFIDENTIAL 


16 


COWPIDENTIAL 


Seotlon  2»     Aasumptlona   leading  to   the  hydraulic  theory. 

We  make   eertaio  simplifying  assumptions   eoncerning 
tba  motion  of  the  water  in  a   tube  of  circular  cross-section. 
The  area  of  the   cross-section  may  vary  along   the   tuhe^   and 


Figure  2.1 

the  vater  is  assumed  to  be   in  motion  ruider  the   pressure   of 
a  volume  of  gas,   as   indicated.      The  following  assumptions  are 
made  i 

1*     The  gas   is  at  every  instant  in  a  uniform  state 
throughout  its  volume,  and   the   interface  between  gas  and  water 
is  a  plane  normal  to  the  axis  of  the   tube*      The  gas   is  assumed 
to  expand  adiabatically. 

2*     The  water  is   considered  to  be  incompressible  and 
non-viscous,   and   its  motion   to  be  irrotational* 

3*     The  hydraulic  assumption  is  made,   i.e.   that  the 
motion  can  be  described  with  sufficient  accuracy  in   terms  of 
the  average  velocity  over  each  cross-section. 

4.     The  pressure  at   the  exit  is   a  constant,   independent 
of  the  exit  velocity  of   the  water. 

The  first  three  assumptions  seem  quite  reasonable,   at 
any  rate  for  the  purpose   of  an  investigation  having  a  prelimi- 
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navy   character.   Various  refinements  could  and  will  be  made 
If  a  need  for  them  should  arise  later.   The  fourth  assumption 
is  a  more  doubtful  one  and  may  prove  to  be  a  rather  poor 
approximation  in  the  cases  of  practical  importance.  We  have 
nevertheless  made  all  of  our  calculations  on  the  basis  of  this 
assumption.  It  seems  reasonable  to  suppose  that  we  can  at 
least  make  valid  comparisons  of  different  cases  since  it  is 
to  be  expected  that  the  inaccuracies  arising  from  the  use  of 
assumption  (4)  affect  the  various  cases  ne   consider  in  the 
same  general  way. 

On  the  basis  of  these  assumptions  it  is  quite  easy  to 
obtain  a  mathematical  formulation  of  the  problem  of 
determining  the  motion  of  the  water,  the  thrust  exerted  on 
the  tube,  and  any  other  relevant  quantities. 
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Section  5>     Fonnulatlon  of   the  hgrdraulic  theory. 

In  this   section  we  derive   the  equation  of  motion  of 
the  water  in  the   tube  baaed  on  the  assumptions  of  the  preceding 
section.      The   Important  fundamental  formulas  for   the   thrust, 
impulse,   and  time   of  stroke  are  also  derived. 


Closed 
End- 


Velocit;^     Ue 


Figure  3,1 

We   shall  assume   that  the   tube  moves  with  constant 
velocity  through  the  surrounding  water.      The  equations  which 
follow  are  written  with  respect  to  a   coordinate   system  moving 
with  the   tube}      the  coordinates  thus  form  an  inertial  system. 

Some  of  the  quantities  which  interest  us  are  indicated 
in  Figure  3.1.     By  Q  and  TJ  we  mean  the   (circular)  cross-section 
area  of  the   tube  and  the  velocity  of  the  water  at  any  section 
of  the   tube.     It  turns  out  to  be   convenient  for  our  purposes 
to  consider  that  any  particular  section  of  the   tube  is  lo- 
cated by  giving  the  volume  V*  of  the   tube  between  that 
section  and  the   closed  front  end  of  the   tube   (instead  of  giving 
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the  distance  between  these  sections,  for  example.)  Thaa  Q 
Is  a  function  of  V*J  Q*  -  Q(V*).  Since  tba  motion  of  the 
water  in  the  tube  is  (in  general)  not  steady,  the  velocity 
U  will  be  a  fimction  of  the  time  t  as  well  as  of  the  volume 
V*:  that  is,  U*  =  U(V*,t). 

The  motion  of  the  water  results  from  the  presstire 
exerted  by  a  layer  of  gaa  located  at  the  closed  end  of  the 
tube.  We  denote  by  V  the  volume  of  the  gaa  at  any  instant, 
by  p  its  pressure,  and  by  V  its  time  rate  of  change.   These 
three  quantities  are  thus  functions  of  the  time  t/   However, 
in  all  of  our  work  we  consider  only  time  intervals  during 
which  the  gas  voliune  is  expanding  (that  is,  motions  such 
that  the  water  is  always  moving  to  the  right  and  out  of  the 
tube)  and  hence  the  time  t  can  be  considered  a  uniquely  defined 
function  of  V.  For  the  derivation  of  our  formulas  it  is 
convenient  to  consider  the  pressure  p  and  rate  of  expansion  V 
as  being  f tactions  of  V.  The  quantity  U*  =  U(V*,t),  the 
velocity  of  the  water  at  any  section,  may  also  be  expressed 
in  terms  of  V  instead  of  t;  ir*=  U(V*,V)« 

We  make  use  of  subscripts  to  denote  special  sections 
of  the  tube.  The  subscript  e  denotes  the  exit  of  the  tube, 
so  that  Pq  and  U^  represent  the  pressure  and  velocity  there, 
while  Vg  is  the  total  volume  of  the  tube  (or  the  value  of  V* 
at  the  ®xit  cross-section,  though  we  dispense  with  the  *  in 
this  case   ).  The  initial  gas  pressure,  gas  volume,  and  rate 
of  change  of  gas  volume  (i.e.  their  values  at  the  time  t  s*  0) 
are  denoted  by  p^^,  V^^,  and  V^^  respectively.  In  all  cases  we 


^  ■ — 

It  should  be  emphasized  that  these  quantities  are  not  fundi  ona 

Later  on  we  are  concerned  only  with  the  volume  V  of  the 
gas  and  can  do  without  the  quantity  V*  entirely.  It  is 
introduced  in  this  section  as  a  necessary  aid  in  deriving 
the  fundamental  equations,  ^ 
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determine  the  motion  up  to  the  time  when  the  gas  has 
expanded  down  to  the  "atmospheric"  pressure  p^.   The  gas 
volume  and  its  time  rate  of  change  at  this  moment  are  denoted 
by  V  and  V  -  the  symbol  p  for  the  gas  pressure  at  this 
instant  is  redundant  since  p„  ^  p^.   The  subscript  s  refers 
to  the  "stroke",  which  is  assumed  to  end  when  the  gas  has 
expanded  down  to  the  pressure  p  • 

We  proceed  to  derive  the  fundamental  equations.  It 
is  convenient  to  base  them  on  the  law  of  conservation  of 
energy.   This  law  expresses  the  equality  of  the  net  work  done 
by  the  external  forces  on  the  original  mass  of  water  in  the 
tube  to  the  change  in  kinetic  energy  of  the  same  mass  of 
water  when  the  gas  in  the  tube  expands  from  the  initial  voltime 
V^  to  some  volume  V.  The  net  work  done  is  then  clearly  given 

ty    (p-P^)dV,  which  represents  the  work  done  by  the  expand- 

ing  gas  diminished  by  tha  work  done  in  forcing  the  water  out 
of  the  tube  against  the  constant  exit  pressure  p  .   (Here  we 
make  our  first  use  of  the  ass\mption  that  the  water  is  incom- 
pressible.) The  change  in  kinetic  energy  is  obtained  by 
subtracting  the  kinetic  energy  K,  of  the  water  in  the  tube 
at  the  initial  instant  (i.e.  when  the  gas  volume  is  V^)  from 
the  sum  of  the  total  energy  flux  K,  out  of  the  exit  section 
(during  the  expansion  from  V.  to  V)  and  the  energy  Kg  of  the 
water  remaining  in  the  tube  at  the  moment  when  the  gas  has 
expanded  to  the  volume  V. 

,      The  basis  for  the  calculation  of  the  kinetic  energies 
is  the  following  formula  for  the  velocity  TJ(V*,V)  at  any 
section  of  area  Q(V  )j 

(3.1)  TJ(V*,V)  ^  V(V)/i(V*), 

which  is  an  immediate   consequence  of  the  assumptions    (2)  and 
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(3)  (see  Section  2  above),  referring  to  the  Incompreasl- 

bility  of  the  water  and  the  use  of  the  average  velocity 

over  any  section.   In  particular,  the  exit  velocity  TJ 

is  given  by  V(V)/^^.  When  the  gaa  volvme  V  changes  by 

an  amount  dV,  a  mass  of  water  pdV  is  forced  out  of  the  exit 

section,  p  being  the  density  of  the  water.   It  follows  that 

the  total  energy  flux  K,  out  of  the  tube  is  given  by 


V 

(3.2)  %  =  "2  ^*  "T  f  ^^^   • 


9  'V^ 


The  kinetic  energy  Kp  of  the  water  remaining  in  the  tube  at 
the  final  instant  is  obtained  by  integrating  wPu  (V  ,V)  with 
respect  to  V*  from  V  to  V^,  V  being  held  fixed:   the  result 
is,  in  view  of  (3.1 )» 


(3.3)  K2-^pv2(V)r    % 


The  kinetic  energy  Kg  of  the  water  initially  in  the  tube  is 
computed  in  the  same  way  as  for  Kg,  with  V  taking  the  special 
value  V.  and  the  limits  on  the  Integral  extending  from 


\   to  V^t 


'^e   « 

®  dV 


(3.4)  S=^^^/   755- 


Q' 


Prom  (3.1)  we  observe  that  the  motion  of  the  water  is 
known  as  soon  as  the  gaa  voliane  V  la  known  as  a  function 
of  the  time  t. 
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The  equation  expressing  the  energy  balance  Is  therefore? 

V 
(3.5)     I        (p-Pg)dV  = 


^ 


1  r    :?..r  ^  ,*r2    r  ®dv*    -2  r  ®dv* 


V  dV   +  V      . 


L^;,        "Jv  ^ "  \  ^ 


■  ''L 


The  differential  equation  of  motion   can  now  he 
obtained  from   (3.5)  by  differentiation  with  respect   to  V. 
In  doing  so,   it  should  be  noted   that  the   second   term  in  the 
brackets  in   (3.5)  furnishes  two   terms  upon  differentiation: 
one  resulting  from  the  factor  v   ,    the   other  from  the  fact 
that  the  lower  limit  of  the  integral  is  V.      The   third  term 
in  the  bracket  is  of  course   independent  of  V.      The  equation 

of  motion  is 

V 

In,  «^     P'^e    .    1  A2/__l             1\      Id  ,*2.    f    ®dV*  ^  o 
(3.6)      — - —   +  w  V    (    tj  n)-  w  ^rr?lV    )•  —m   -  0, 

P  ^  Q    (V)        Q^^      ^   ^^  J^     Ql'^ 

in  which  the  quantity  Q(V)   in  the  parenthesis  refers   to  the 

cross-section  area  at  the  gas-water  interface.      The  equation 

(3.6)  can  also  be  derived  from  the  Bernoulli  law.      The 

differential  equation   (3.6)   is  a   first  order  differential 

equation  for  the  quantity  V,   that  is,  for  the  rate  of 

expansion  of  the  gas.     If  the  time  were  to  be  reintroduced 

by  the  relation  dV  =  Vdt,   the  result  Is,   of  course,  a   second 

order  differential  equation  for  V(t).     In  fact,  w  ;nr(V   )  - 
1    d/;,2.  AV  ^U  ^  ^^ 


So  far  we  have  not  specified  the  pressure  p  of  the  gas. 
Since   the  gas   is  assumed   to  expand  adiabatically  (assumption 
(1))  we   can  express  p  as  a  function  of  V  by  the  relation 

(3.7)  pv''=Pj^Vj^, 
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in  which  y  is  the  adiabatlc  exponent. 

The  equations  (3.6)  and  (3.7)  together  with  the  values 
V^  and  V^  for  the  gas  volvaae  and  Its  rate  of  change  at  the 
initial  instant  constitute  a  complete  mathematical  formulation 
of  the  problem  of  determining  V  as  a  function  of  t  once  the 
shape  of  the  tube  and  values  for  the  exit  pressure  p  and  the 
adiabatic  constant  y   have  been  given.   The  succeeding  sections 
of  this  report  are  concerned  with  the  solution  of  this  initial 
value  problem  and  its  interpretation  for  practical  purposes. 
In  the  remainder  of  this  section  we  obtain  formulas  for 
various  important  quantities,  which  can  be  computed  once  the 
differential  equation  has  been  solved  for  V  as  a  function  of  t* 

Two  of  the  most  important  basic  quantities  of  interest 
on  practical  grounds  are  the  impulse  and  the  thrust  made 
available  during  the  stroke  for  propulsive  purposes.  In  the 
above  derivations  we  have  fixed  our  attention  on  the  forces 
exerted  on  the  water  in  the  tube  rather  than  on  the  forces 
exerted  on  the  tube  itself.  Since  we  wish  to  continue  our 
analysis  from  this  point  of  view  because  of  its  convenience, 
it  is  necessary  to  indicate  the  relation  between  the  net  force 
and  the  net  impulse  imparted  to  the  water  and  the  net  forward 
thrust  and  Impulse  imparted  to  the  tube.  In  all  that  follows 
the  quantities  F  and  I  represent  the  resultant  force  on  the 
water  in  the  tube  at  any  time  and  the  total  impulse  delivered 
to  the  water  in  the  tube  from  the  initial  Instant  up  to  any 
time.  Both  of  these  quantities  are  positive  when  directed 
"backwards",  i.e.  in  the  direction  toward  the  exit  of  the  tube* 

The  following  Figure  3.2  indicates  the  relation 
between  F  and  the  forces  on  the  tube. 
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(a)  P  =  total  forces  on  water 


("b)  D  =  total  drag 

Figure  3.2 

The  force  F  on  the  water  arises  from  the  forces  exerted  by 
the  tube,  the  gas,  and  the  water  pressure  p^  at  the  exit,  the 
latter  being  directed  opposite  to  P.  Consequently  the  tube 
exerts  a  backward  force  P  +  p_Q-  on  the  water  contained  In  It 
The  total  force  on  the  tube  arising  from  the  water  Inside  it 
la  thus  directed  forward  and  has  the  magnitude  F  +  p^Q^.  If 
the  average  pressure  of  the  water  on  the  outside  of  the  tube 
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is  greater  than  p  ,  the  excess  of  the  average  pressure  above 
p  will  create  a  drag  which  we  include  in  the  total  drag  D 
(see  Figure  3.2b)  along  with  the  other  components  of  the 
resistance  arising  from  the  hydropulse  and  the  object  it 
propels.  Consequently  we  see  that  the  quantity  F  we  have 
defined  represents  numerically  the  thrust  available  to  over- 
come the  drag.  In  the  same  way  the  quantity  I  (which  is  of 
course  the  time  integral  of  P)  is  seen  to  represent  the 
net  impulse  available  for  propulsive  purposes. 

We  obtain  first  a  formula  for  the  impulse  I,  regarded 
as  a  function  of  the  gas  volume  V.  By  the  law  of  conservation 
of  momentum,  it  is  equal  to  the  momentum  of  the  water  remaining 
in  the  tube  at  the  final  instant,  plus  the  momentum  trans- 
ported through  the  exit,  minus  the  momentum  in  the  water  at 
the  initial  instant.  By  the  same  considerations  as  were  used 
in  deriving  (3.5)  we  find 

(3.8)  I(V)  =  p  \^(   VdV  .   V  f '^  -  KI"'—^ 

It  has  already  been  implied   that  the   time   interval  t 

for  the   complete  exhaust  stroke  is   chosen  as   that  extending 

from  the   initial  instant  t  =  0  to   the   Instant  when  the  gas 

voliome   is  V     (or  the  gas  pressure   is  p..)      The   total  impulse 

per  stroke  I   ,  which  is  a  fundamentally  important  quantity 
s 

for  us,   is  obtained  from  (3.8)  by  setting  V  -  Vg.     That  lat 

(3.9)  Ig  =  I(Vg). 


4 

Since  we  assume  the  tube  to  be  moving  with  constant 

velocity  we  have  P  =  D. 
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The    thrust  P  can  be  obtained  by  differentiating  the 

dV       • 
impulse   I  with  respect  to  the   time.     We  observe   that  gr  =  V 

2  •      • 

and  that  ^-^  =  V  ^§|=V^and  obtain  from  (3.8)  the 


following  formula  for  the  thrust: 

(3.10)      P  =/>v2(^  -  ^)  ^    pV  §^J 


^e  * 
®dV 

■J'  • 


e   '*^'"       "'-^y   Q 

It  is  of  great  interest  to  compare  the  relative 
quantities  of  fuel  required  for  the  exhaust  stroke  under 
various  operating  conditions.  This  should  be  done  by 
considering  in  detail  the  reaction  of  the  fuel  with  the  water, 
but  since  we  have  little  or  no  information  on  this  point  we 
have  made  our  calculations  referring  to  fuel  consumption 
assuming  that  initially  we  have  in  each  case  a  layer  of  an 
Ideal  gas  of  a  specified  molecular  weight  at  a  definite 
temperature.  Thus  we  fix  the  mass  M  of  fuel  used  in  any 
given  case  by  the  formula* 

(3.U)  Mg  ^  5^  p,V^  . 

in  which  fi  ia   the  molecular  weight,  R  is   the  gas  constant  per 
mol,  6t    is  the   initial  absolute   temperature,  and  p.   and  V 
are   the   initial  pressure  and  volume  of  the  gas   layer. 

One  of  our  principal  objects  la   to  study  the  dependence 
of  the  impulse  per  stroke  per  unit  mass  of  fuel,   that  ia,   the 
quantity  Ig/^e'  °°  variations   in  all  of  the  relevant  parameters. 
The  quantity  Ig/^a  ^^   ^^®  dimensions  of  velocity  and  is 
called  the  effective  velocity  (or  specific  impulse.)     If  we 

set  Ig/M     =  H~7T  ^^^li  ^  the   time  for  a  complete   cycle,   then 
we   see   that  the  effective  velocity  represents  the  ratio  of 


CONFIDENTIAL 


27 

CONFIDENTIAL 


average  thrust  to  rate  of  fuel  consumption  (provided  we  neglect 
the  losses  in  impulse  during  the  refill  and  scavenging  portions 
of  the  full  cycle.)  Thus  the  effective  velocity  is  a 
reasonable  measure  of  the  efficiency  of  the  hydropulae.  We 
shall  study  its  dependence  on  variations  in  initial  gas 
pressure,  mass  of  fuel,  initial  velocity  of  the  water  in  the 
tube,  and  on  variations  in  the  adiabatic  exponent  y  •  We 
shall  also  study  the  dependence  of  the  time  of  stroke  t  , 

and  the  average  thrust  over  the  stroke  "F  *=  I  /t  on 

'  ■  '       ^ '       '  '  s  3/3 

variations  in   the   same  parameters. 

The  time  of   stroke   t     is  obtained  from  the  formula 

3 


(3.12)  t^  =J  °  3V  , 

in  which  V  is  to  be  determined  by  integration  of  the 
equation  of  motion  (5«6)« 
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Section  4>     Dlmenalonlesa  variables  and  parameters. 

For  certain  of  the  later  sections   It  Is   convenient  to 
make  use   of  dlmenslonless  variables  and  parameters,   particular- 
ly «hen   the  actual  Integration  of  the  equation  of  motion   (3.6  ) 
Is   In  question. 

For  the  dlmenslonless   time  r  we  write 

(4.1)  t  =  T  •  T    , 

in  which  the  unit  of  time  T  Is  defined  by 

(4.2)  T  ;=-  ® 


All  quantities  In  (4.2)  have  been  defined  In  the  preceding 
section.  The  dlmenslonless  gas  volume  v  Is  defined  by 

(4.3)  V  =  V^   •  V   , 

t 

while  the  dlmenslonless  rate  of  expansion  v     la  given  by  the 


relation 


(4.4) 


V  =  Qe/-X  ^  • 


The   prime  denotes  differentiation  with  respect  to   r  •     We 

observe   that  the  quantity ^—2^  has   the  dimensions  of  velocity. 
It  la  perhaps  of  interest  to  note   that  the  value  of  this  unit 
of  velocity  is  about  32  ml./hr.  when  p     is  atmospheric  pressure. 

*  In  certain  later  sections  a  different  unit  of  velocity  la 
sometimes  used,  but  this   occurs  only  in  formulas  not  written 
in  terms  of  the  dlmenslonless  quantities  defined  here. 
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We  denote  by  q  the  dlraenaionlesa   croaa-aection  area  defined  ty 

(4.5)  Q  =  Qq   •  q  . 

Finally,  we   Introduce   the  dlraenslonlesa  quantltlea 
u,  f ,  and  j   corresponding  to  the  velocity  of  the  water  U, 
the  thrust  F,  and   the   impulse' H 

(4.6)  U  =  /— 2.  .  u  , 

V    P 

(4.7)  p  =  2p^Q^    .   f   , 

(4.8)  I  =PV^J-^  .    j   . 

In  terras  of  our  dlmenalonless  quantities  the  differential 
equation  of  motion   (3.6)  derived  In  the  preceding  section 
takes  the  form 

1 

(4.9)  vy   .  1  ^  (v*)2(J^  -  1)   -  ^(v*)2    \    %^0, 

^  V   q 

In  which  Vg  =  Vg/V^  and  the  prime  refers  to  differentiation 
with  respect  to  r  • 
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Section   5.     Duct  of  variable   croaa-aectlon.     Approximate 
formulas.     Shape  of  duct  for  maxlmvun  Impulse. 

The   Integration  of  the  differential  equation  of  motion 
(3.6)   Is  not  difficult,  but  very  laborious,   since  it  must  be 
done  numerically.      Consequently  a  complete  discussion  of   the 
solutions  of  the  differential  equation  and  of  the  important 
quantities,    such  as   the   total  Impulse   and  the   time   of  stroke, 
in'  their  dependence  on  all  of  the  many  parameters  would  involve 
a  mountainous  labor  in  calculation  if  it  were  necessary  to 
carry  out  numerical  integrations   in  every  case.     Fortunately, 
it  turns  out  that  all  of  the  more   important  quantities  can  be 
given  by  relatively  simple  and  sufficiently  accurate  approxi- 
mate formulas,  which  contain  all  of  the  parameters  explicitily 
and  thus   lend   themselves  well  to  a   complete  discussion. 
These  formulas  have  been  checked  by  comparison  with  the 
results  obtained  by  numerical  integration  of  the  differential 
equation  in  quite  a  large  number  of  cases   (including   two   types 
of  tubes  with  conical  front  portions,   as  well  as   straight  tubes) 
so  that  limits  for  the  accuracy  of  the  approximate  formulas 
are  known  (and  will  be  given  later.)     We  proceed  to  derive 
the   approximate  formulas. 

The  basis  of  the  approximate  formulas  is   the  fact, 
observed  in  many  cases  which  were   computed  exactly,   that  the 
pressure  in  the  gas  decreases  very  rapidly  initially  and   then 
stays  relatively  constant.       At  the  same  time,   the  mass  of 
water  in  the   tube   is  strongly  accelerated  initially  and  quickly 
reaches  a  high  velocity  which  then  remains  nearly  constant 
over  the  remainder  of  the  stroke.     This  indicates   that  the 
motion  of  the  water  may  be  analyzed  with  good  accuracy  on  an 

This  results  largely  because   the  gas   is  assumed  to  obey 
the  adiabatic  law  so   that  the  pressure  decreases  rapidly 
with  increase  in  volume. 
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impulse  basis  simply  by  assuming  that  the  mass  of  water 

involved  remains  constantly  equal  to  the  initial  mass  in  the 

tube. 

We  denote  by    £     the  energy  imparted   to  the  water  as 

s 

the  net  result  of  the  work  done  by  the  gas  pressure  and  the 
pressure  at  the  exit  of  the  tube;  thus  £  ia  given  by 


/ 


^a 


(5.1)  ^a  "  J    ^P  -  Pe^  ^^  • 

in  which  V.  and  Vg  represent  the  initial  and  final  volumes 
of  the  gas,  V_  being  the  gas  volume  at  the  end  of  the  stroke 
when  the  gas  pressure  has  the  value  p  .  We  denote  by  M  the 
initial  mass  of  water  in  the  tube  and  by  £,  and  Af.    the 
kinetic  energy  and  momentum  of  the  water  in  the  tube  at  the 
initial  instant.  We  then  use  the  energy  law  to  express  £ 

0 

by  the  following  approximate  formula: 

1    -2 

(5.2)  ^3  -S^w  ^a  -  ^1  ' 


in  which  U_  is  the  average  of  the  square  of  the  final 
s 

velocities  of  the  water  particles  in  the  tube.  The  impulse 
per  stroke  I  —  a  fundamentally  important  quantity  for  us  — 
can  be  approximated  as  follows} 

(5.3)  I3  =M,  U3  -  Af^ 


It  should  be  clear  that  our  approximation  refers  to  the 
quantity  I3  as  it  was  defined  in  Section  3. 
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In  which  U     is   the  average  of  the  final  velocity  of  the  water 
s 

particles  in  the   tube.     We  now  make   the  assumption  that 

2 

U  =  U^,  I.e.  that  the  average  of  the  squares  of  the  velocities 

s     8' 
Is  equal  to  the  square  of  the  average  velocities,  (this 

assvunptlon  is  exact  for  the  case  of  a  tube  of  constant  cross- 
section)  and  eliminate  TJ„  and  U„  from  (5.2)  and  (5.3)  to 

s  s 


obtain  the  important  approximate  formula 


(5.4)  Ig  «     /2M,£g   +  2M/j_     -    Af^ 

for  the  total  impulse  per  stroke   in  terms  of  quantities  which 
can  be   calculated  or  estimated  in  any  given  case  without  much 
difficulty.     An  important  special  case   is   that  in  which  the 
initial  velocity  TJ^^  of  the  water  in  the  tube  is  zero,  and  I^ 
is  given  by  the  simple  formula 


(5.5)  Ig  =  /2M,£. 


w 


s  • 

The  formula  (5.4)  can  be  Interpreted  as  the  result  of 

replacing  the  hydropulse  by  a  simpler  mechanical  modelJ   The 

formula  furnishes  the  total  Impulse  on  a  rigid  body  of  mass  M. 

which  moves  in  a  straight  line  under  the  action  of  forces 

which  do  work  of  amount  £_.   The  quantities  TJ_  and  U  in  (5.2) 

s  s      s 

and  (5.3)  are  then  both  equal  to  the  square  of  the  final 
velocity  of  the  rigid  body. 

The  time  of  stroke  t_  Is  determined  by  the  approximate 

s 

formula 

(5.6)  tg^iliZji^, 

s 
In  which  Q_   is  the   cross-section  area  of   the   tube  at  the  exit. 
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and  Ug  is  to  be  determined  from  (5.3).   That  (5.6)  furnishes 

some  sort  of  approximation  to  the  correct  value  for  t  is 

s 
clear. 

an 

Finally,    the   average    thrust  P_   on   the  water  in   the    tube 

3 

during  the  stroke  is  given  approximately  by 


(5.7) 


^3  =  Is/* 


3/  3 


With  Ig  and  tg  to  be  determined  by  (5.4)  and  (5.6). 

It  is  important  to  know  that  the  approximate  value 
given  by  (5.4)  for  Ig  is  a  rigorous  upper  bound  for  this 
quantity.  We  proceed  to  prove  this  statement.  Prom  the 
formulas  (3.8)  and  (3.5)  of  Section  3  we  have 


(5.8) 


and 


I3  ^M^   =p 


VdV  +  V. 


^3    - 


(5.9) 


E     +  E 

^3    ^i 


h 


1 

Lre 


•2 

V^  dV  + 


^3/ 


«dV* 


We  shall  write  the  two  integrals  in  (5.8)  as  one  single 
integral,  first  writing  the  second  of  the  integrals  in  (5.8) 


as 


®dV 
-^  ,   which  is  legitimate  since  the  variable  of  inte- 


£1 

gration  is  a  dummy  variable.  Since  the  intervals  of  Integration 
in  the  two  integrals  (3.3)  are  contiguous  and  do  not  overlap, 
we  may  write 

(5.10) 

a  i.  J 


^±  "^  '^J   H(V)dV 
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where 


r      * 


H(V)   = 


V/Qq     for     Vj^   <  V  <  V^ 


L   Vg/^     for     Vg   <  V  <  V^ 


Similarly, 
(5.11) 


£g   +  E^  -  ^pj      H(V)^dV  . 
^i 


Applying   the   Schwarz   inequality  to    (5.10),   we   have 
(5.12)        Ig   +    /J/^   <   P  (Vg   -  V^)V2     I       H(V)^dV 

Lv,      J 

(The  Schwarz  inequality  asserts  that 


I 


f(V)g(V)dY  </  I  f(V)^dV  •      g(V)^dV   , 


equality  holding   only  if   the  ratio  of  f(V)   to  g(V)   is   a 
constant.      To   apply  it   to    (5.10),    take  f(V)   =  1     g(V)   -  H(V)    ). 

Using    (5.11)  and   the  fact  that  the   initial  mass  M^  of 
the  water  in   the   tube   is  M_  -  P(V^  -  V.),  we   obtain 


^8  -^  ^i  :^   ysM^Cf^g  +£^i) 


or 


( 5. 15 ) 


Is^    /2M^£g   +2M,f:^'-  Mi, 


In  particular,  if  the  initial  velocity  of  the  water  In  the 
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tube   la   zero,  we  have 

(5.14)  I     <    v/2M„  E^    . 

Thus,  the  approximate  formulas  (5.4)  and  (5.5)  yield  rigorous 
upper  bounds  for  I  • 

It  is  now  possible  to  draw  a  number  of  conclusions  of 
some  practical  Importance  regarding  the  effect  of  changes  in 
shape  of  the  hydro pulse  tube  on  the  Impulse  per  stroke.  For 
this  purpose  It  Is  convenient  to  write  the  inequality  (5.13) 
In  the  form 

(5.13)'        Ig  +  ^i  <  /2M^  Eg  +  2M^  E^   , 

and  to  recall  (cf.  equation  (5.10))  that  Ig  +Af.    =P    H(V)dV, 

and  that  (5.13)  results  from  application  of  the  Schwarz 
inequality  to  the  integral.   Our  calculations  for  many  cases 
show  that  the  upper  bound  for  Ig  +  A/j,  furnished  by  (5.13)  is 
practically  attained  for  the  straight  cylindrical  tube  (cf. 
the  table  at  the  end  of  the  following  Section  6)  as  well  as 
for  a  number  of  tubes  having  conical  portions  (see  Section  8). 
The  reason  for  this  lies  in  the  behavior  of  ^  as  a  function  of 

V  over  the  stroke:   Typical  cases  for  the  straight  tube  are 
shown  graphically  in  Figure  7.9  of  Section  7,  which  shows  that 

V  increases  very  rapidly  at  first  but  then  remains  constant  — 
in  other  words,  V  is  practically  constant  over  most  of  the 
range.  Hence  the  function  H(V)  also  behaves  like  ^a  constant 


In  Sections  6  and  8  the  dimensionless  variables  defined 
In  Section  4  are  used  rather  than  the  variables  of  the 
present  section. 
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(cf.  the  definition  following  equation  (5.10))  if  the  cross- 
section  area  of  the  tube  does  not  vary  too  much  so  that 
the  Schwarz  inequality  is  practically  an  equality. 

Let  us  consider  first  cases  in  which  the  initial 
velocity  of  the  water  in  the  tube  is  zero,  so  that  the  initial 
momentun  Mm   Is  zero.  In  this  case  formula  (5.14)  applies, 
of  course. 

It  can  now  be  seen  that  variations  in  the  shape  of  the 
tube  in  these  cases  have  only  a  slight  influence  on  the  impulse 
per  stroke,  at  least  for  tubes  which  are  not  too  radically 
different  from  the  straight  tube:   This  is  an  immediate  conse- 
quence of  the  fact  that  the  upper  bound  for  I  given  by  (5.14) 
is  practically  attained  for  the  straight  tube,  while  on  the 
other  hand  the  right  side  of  (5.14)  is  independent  of  the 
shape  of  the  tube  —  it  depends  only  on  the  volume  occupied 
by  the  original  mass  of  water  M„  in  the  tube  and  on  the  energy 
imparted  to  this  mass  (which  depends  only  on  the  initial 
pressure,  volvune,  and  character  of  the  gas  and  not  on  the 
shape  of  the  tube.)  Thus  if  M_  is  kept  fixed  the  straight 
tube  will  yield  very  closely  the  maximum  impulse  with  respect 
to  variations  in  shape.   The  determination  of  the  exact  shape 
for  maximizing  the  impulse  is  therefore  settled  in  case  the 
forward  velocity  of  the  hydropulse  is  not  too  great. 

The  cases  in  which  the  initial  velocity  of  the  water 
in  the  tube  is  not  zero  are  not  so  easily  dealt  with.   Of 
course,  if  the  Initial  velocity  is  not  too  great  (i.e.  if  the 
forward  velocity  of  the  hydropulse  is  not  too  great)  then  our 
conclusion  that  the  straight  tube  is  as  good  as  any  other  for 
the  purpose  of  maximizing  the  impulse  would  still  be  reason- 
able. For  high  velocities,  however,  the  quantity  M.    in  (5.13) 
appears  to  be  sensitive  to  changes  in  shape  of  the  tube.  In 
Section  8,  where  a  few  cases  of  tubes  with  conical  portions 
are  treated,  the  Impulse  for  an  initial  velocity  of  about  90 
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ml.  AiT   is   increased  by  from  10  per  cent  to  20  per  cent  by 
taking  a   conical  tube   in  place  of  a  straight  cylindrical   tube. 

For  the  higher  velocities  of   the  hydropulse  we  have 
not  settled   the  question  of   the  best  shape  of  tube  for  maxi- 
mizing  the   impulse,   and  more  work  should  be  done   on   this 
problem.      However,  we  have   obtained  one  further  partial  result 
which  is  not  without  interest:     Consider  only  that  part  of 

the   tube   between  V„  and  V^,    i.e.    that  part  which  is  still 

s  a ' 

filled  with  water  at  the  end  of  the   stroke.     By  applying  the 

Schwarz  inequality  to   the   second  integral  on  the  right-hand 

side  of   (5.8)   it  can  easily  be   shown  that  the  impulse  is  a 

maximum  when   this   part  of  the   tube  has   a  constant  cross-section, 

if  the  volume   of  this  part  is  assumed   to  be  fixed  in  all  cases 

under  comparison.        This  result,   together  with  the  fact  that 

the  calculations  of  Section  8  show  the   impulse  at  high  initial 

velocities   to  be  increased  by  flaring  out  the  front  section, 

indicates   that  the  best  shape  of  tube   for  maximizing   the 

impulse   is  one  with  a  flared  out  front  section  attached   to  a 

straight  cylindrical  portion  extending   to   the  exit  section. 

We  proceed  next  to  show  that  th©  approximate  formula 

(5.6)  for  t     furnishes  a  rigorous  lower  bound  for  this 

quantity  in  the  case  of  th©   tube  of  constant  cross-section. 

It  also  undoubtedly  furnishes  a  lower  bound  for  tubes  which 

do  not  depart  all  too  radically  from  the  straight  cylindrical 

tube,  but  we  do  not  discuss   this  question  here.      (See,  however. 

Section  8  where   t       is  calculated  for  two  different  types  of 

conical  tubes,  which  depart  rather  radically  from  the   straight 

shape.      The  approxifliate  formula  yields   t     too  low  in  these 

s 

Upon  what  basis  the  comparison  should  be  made  is  not 
entirely  clear.  This  point  Is  discussed  in  Section  8. 

Without  such  a  restriction  the  total  impulse  could  easily 
be  seen  to  have  no  finite  upper  bound. 
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cases  also.)  We  begin  with  the  exact  formula  (3.12)  for  t  : 


(5.15) 


J..   V 


We  write   (5.15)   in  the  following  form  by  adding  a  certain 
integral  to  both  sides; 


Vg  -^v^     e  e 


1  s 


in  which  Q  =  Q(V  ),  as  defined  in  Section  3.   Through 

introduction  of  the  ftmction  H(V)  defined  above  (cf.  equation 

(5.10)  and  the  remarks  which  immediately  precede  it)  and  a 

division  of  both  sides  of  (5.16)  by  V^-  V.  we  obtain  the 

e   1 

relation 


Q 


e  1 


*a* 


^1 

V  J 


e 


V  J   **e 


dV 


'"^t 


V. 


in  which  the  dummy  variable  V  has  been  replaced  by  V.  A 
known  inequality  asserts  (see  the  book  of  Hardy,  Llttlewood, 

Polya  on  Inequalities,  p.  151)  that  if  0(y)  is  convex  — 

w 

I.e.   if  0     (y)  >  0  —  then 


a 


))dV  >    0 


^    J 


f(V)dV 


Applying   this   inequality  to   the  right-hand   side   of   (5.17) 
with  0(y)   =i  and  f(V)  =  H(V),  we  obtains 

ft 
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(5.18) 


Q 


e 


*s  * 


3 


ti 


e 


Q 


dV     >  - 

J  7. 


o 


H(V)dV 


By  making  use  of  (5.10)  (an  exact  formula,  it  should  be  noted) 

to  replace  I   H(V)dV  and  noting  that  p(Vq-V^)  =  M^  we  obtain 

^i 
from  (5.18)  the  result 

s 

We  know  from  our  calculations  that   the   lower  bound 

furnished  by  the  right-hand  aide   of  (5.19)   is  practically 

attained  for  all  tubes  not  too  radically  different  from  the 

straight  tube.     Unfortunately,   the  right-hand  aide  of   (5.19) 

contains  the  quantity  V_  which  can  be  determined  only  by 

s 

integration  of  the  equation  of  motion.  However,  in  the  most 
important  case  —  that  of  the  tube  of  constant  cross -section  -  ■ 
it  is  possible  to  eliminate  V  ,  as  follows; 

We  know  that  V  is  an  increasing  function  of  V  for 
V^  <  V  <  V  .  Prom  (5.10)  we  may  therefore  write  in  this  case 
(Q(V)=Qg): 


(5,20) 


Ig  +  '^l  <  P  V  .  -V^  ,  or 


^ 


e 


(5.21) 


V^-V. 


-L  £.       _e  '^i 
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Substituting  from  (5.21)  In   (5, 19)  Tie  obtain 

or,   since    p  (V^-Vj^)  =  M^: 


8 


Since  (I„  +  Mj,)AS.^   =  TJ„  (from  (5.3)),  i»e  see  that  the  right- 

S       X  '  w     3 

hand  side  of  (5.23)  Is  the  same  as  that  of  (5.6).  Hence  (5.6) 

furnishes  a  rigorous  lower  bound  for  t  in  the  case  of  the 

s 

tube  with  tiniform  cross -section. 

If  I_  and  t     are   calculated  using   (5.4)  and   (5.6)  we 
s  s 

may  now  conclude   that  the  formula   (5.7)  for   the  average   thrto  t 

P     yields  a  rigorous  upper  bound  for  the  straight  tube,   since 
s 

t     is  a  lover  boxrnd  in   this  case  and  I^  is  in  all  cases  an 
s  s 

upper  bound. 

The  acctiracy  of  the  approximate  formula  (5.6)  —  i.e. 
of  (5.23)  using  the  sign  of  equality  ~  is  discussed  at  the 
end  of  th9  following  Section  6.  In  this  case  the  formula 
is  found  to  be  more  accurate,  that  is,  the  lower  bound  la  more 
nearly  attained,  when  the  initial  velocity  is  high  rather 
than  low  (the  opposite  was  true  for  the  impulse).  The 
approximation  is  good  enough  for  practical  purposes. 
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Section  6»      The   tube  of  constant  cross-section*     Formulas 
showing  dependence  of  essential  quantities  on  all  parameters* 

In  this  section  the  approximate  formulas  of  the 
preceding  section  are  adapted  to  the   special  case  of  the   tube 
of  constant  cross-section.      The   integrations  are   carried  out 
and   the  formulas  are  written  in  such  a  way  as   to  make   the 
discussion  of   the  dependence   on  the  various   parameters   a 
relatively  simple  matter.      The  actual  discussion  is   then 
carried  out  in  the  following  Section  7.     The  essential  quanti- 
ties for  which  formulas  are  obtained  are:      1)      the  effective 
velocity  Ig/*I-     2)     the   time  of  stroke   tg,   and     3)     the 

average   thrust  P„  over  the   stroke.     In  each  of  these  formulas 

s 

the  dependence   on  all  of  the  following  parameters  is  exhibited: 

1)     the   "fuel  parameter"  M   /M^,   i.e.    the  ratio  of  the  mass  of 

gas  per  stroke   to   the  mass  of  water  in   the   tube,     2)     the 

initial  gas  pressure  ratio  V>^/Pq>     3)      the   initial  velocity* 

U,    of  the  water  in  the   tube,     4)      the  adiabatic  exponent     y« 

We  begin  by  calculating  the  quantity    E.  in  the 

s 

formula  {5.4)«  Por  It  wo  have  (from  (5.1)]: 

(6.1)  £.  =  J   (-Li  -  pj  dV  , 

from  the  adiabatic  gas  law.   The  integration  can  be  carried 
out  to  yield 

(6.2)  £,  =  p,V,   [^(1  -  (^)>'-l)   -  ^i^  .  in    . 

S  X       1 


« 


This  velocity  might  be   considered   to  be  approximately 
equal  to  the  forward  velocity  of  the  hydro pulse. 
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Again  using  the  adiabatic  gas  law  we  may  write  for  Vg  the 
equation 

(6.3)  Pe^s^Pi^i- 

Elimination  of  V„  from  (6.2)  yields 

(6.4)        £3  =  P.v.c^iTd '—r  ^  vjp;'^  • 

(Pi/Pe)  y 

The  mass  M  of  the  gas  is  related  to  p^^V^  by  the  ideal  gaa 
law  (3.11);  replacing  V^^^   ^  means  of  this  relation  we 
obtain  finally  for  £„  the  formula 

(PlA'e^   V 

The  formula  (5.4)  for  the  impulse  per  stroke  I„  can 
be  written  in  the  form 


(6.6)        I3  =  (13)^ 


with 


(6.7)  ^^s^o^sj^^' 

The  quantity   (!_)«  represents   the  Impulse   per  stroke   in  the 

3  O 

special  case  in  which  the  initial  velocity  of  the  water  in 
the  tube  is  zero.  For  the  effective  velocity  Ig/^-  w©  may 
therefore  write 
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(6.8) 


with 


(6.9)      (Ig) 


!./"«   = 


-ifrr 


Pe         ^         P9 


If  Tue  were    to  divide  equations    (6.8)  and    (6.9)  on 

both  sides  byi/S— i  ,  which  has  the  dimensions  of  a  velocity, 
they  would  be   in  dimensionless  form.     All  the  essential 
parameters  are  exhibited  explicity  in  these   formulas. 

It  is  convenient  to  introduce  si)ecial  symbols  for  two 
combinations  of  quantities  in  the  formula  (6.8),  as  follows; 
We   set 


(6.10) 


Pe         ^         Pe 


)+-£-,   and 
Pi 


(6.11) 


and  obtain  the  following  relation  for  the  effective  velocity, 
replacing  (6.8): 


(6.12) 


h;=\ 


R0J 


i-Si-j/^^-kj. 


The  advantage  of  this  form  is  that  h  contains  (aside  from  7  ) 
only  the  pressure  ratio  P^/P-*  while  k  is  the  only  quantity 
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lihlch  involves  the  initial  velocity  TJj^. 

We  proceed  next  to  amplify  the  formula  (5.6)  for  the 

(V„-V,  )>Q^ 
time  of  stroke  tgJ  tg  =   -■  y-  •  •  -  for  the  case  of  the  tube 

of  uniform  croas-aoction.   Thii  is  done  hy  following  much  the 

same  procedure  as  above,  including  the  introduction  of  the 

quantities  h  and  k  defined  by  (6.10)  and  (6.11).   The  result 

for  t  is  readily  found  to  be: 


s 


(6.13)       t,  =   ®'  ^ 


\^  dY' 


3 


W 


/ 


Finally,  the  average  thrust  Pg  over  the  stroke  is 
given  by  the  quotient  of  Ig  and  tg.   The  result  is 

Pi 
(6.14)  P,  =  ^  =  2p  Q  .  -JS—  .  J^^iJ^^   -  k) 


-^s         '-^^        ..  1 


(sr-' 


The  formulas  (6.10)  to  (6.14)  inclusive  make  an  almost 
complete  discussion  of  the  ideal  performance  characteristics 
of  the  hydropulse  possible.  They  do  not,  of  course,  make  it 
possible  to  discuss  the  motion  of  the  water  as  a  function  of 
time.  Also,  these  formulas  are  approximations,  though  rather 
accurate  ones.  In  fact,  the  effective  velocity  is  in  error 
by  less  than  3  per  cent  if  the  pressure  ratio  Pjl/Pq  is  greater 
than  10.  The  time  of  stroke  t„  and  the  average  thrust  P_  are 
subject  to  an  error  of  about  5  i)er  cent  if  Pj^/Pg  >  10  and  the 
initial  velocity  U^  is  greater  than  about  30  miles  per  hour; 
if  TJ^^O  the  error  is  not  greater  than  15  i)er  cent. 

The  curves  of  Figure  6.1  indicate  the  accuracy  of  the 
approximate  formula  for  the  impulse  in  the  case  y  =   1. 5 ,  U.  =  0. 
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The    curves  give    the  values  of  the   dlmensionless    Impulae    j 
(defined  by  equation   (4,8)   of  Section  4  for  various   values 
of   the   initial  pressure  ratio  P^/Pq  as  functions  of   the 
initial  volume   v,    =  V./V  •      The   points   indicated  by   circles 
were   obtained  by  numerical  integration  of   the  equation  of 
motion  and  are   subject  to  an  error  of  not  more   than  1  per  cent. 
The   curves,   on   the   other  hand,   were    computed  by  using    {6»8). 
As  we   see,    the  approximate   formula   furnishes  quite   accurate 
results  over  a  wide  range  of  initial   pressures   and  initial 
gas  voliaraes.      The  following   table   gives  a  number  of  additional 
comparisons  for  j     when    7    and  U.    are  varied.      (The  quantity 

,31 

v^    is   the  dimensionless  initial  velocity  defined  in  Section  4, 

••■  t 

If  p_    is  atmospheric  pressure,   then  v^   =  1  corresponds   to  a 

velocity  of  31.8  miles/hr.) 


Accuracy  of  Approximate  Formula 
for  Impulse 


Cases 

Impulse  i 
s 

Approximate 

Exact 

^  =  ^•5   V,  =.01  'i=l 
v^  =  .05,  Vjj^  =  3 

1.56 
.82 
.069 

1.57 
.84 
.068 

7  =  1.1            '1  ~  ° 

V  =  0.5 

°                                           f 
v^=  .01           V.  =  3 

1.373 
.302 

1.366 
.301 

7  =  1.3            v^  =  0 

v„=  .416 

v^=     .024          ^1  ~  -^ 

1.198 
.570 

1.188 
.566 
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A3  the  table  shows,  the  approximate  formula  appears  to  be 
quite  accurate  over  a  wide  range  of  parameter  values. 

In  the  following  table  a  few  conparlsons  of  exact 
and  approximate  values  for  the  time  of  stroke 
dimensionless  form,  cf.  Section  4)  are  given: 


T3  (in 


Cases 

Time  of 

Stroke  Tg 

Approximate 

Exact 

7  =  1.1  Vj^  =  .01 
Vg  =  0.5 

v;.o 
v;  =  3 

.353 
.148 

.392 
.150 

7  =  1.3  v^   =   .024 
V3  =  .416 

v;  =  o 

vl=l 

.319 
.248 

.370 
•  259 

As  remarked  above,  the  approximate  formula  for  the  time  of 

t 
stroke  becomes  more  accurate  when  the  initial  velocity  (v. 

in  dimensionless  form)  increases. 

As  we  have  seen  in  the  preceding  section,  our  formula 

for  I_  gives  an  upper  bound  for  this  quantity.   The  formula 
s 

for  t_,  as  we  have  seen,  furnishes  a  lower  bound,  and 
consequently  the  formula  for  P_  yields  an  upper  bound. 

The  following  section  is  devoted  to  a  discussion  of  the 
formulas  developed  here,  largely  on  the  basis  of  graphs.  We 
also  include  some  results  obtained  by  numerical  integration 
of  the  differential  equation  —  in  particular,  the  gas  pressure 
is  given  as  a  function  of  time  in  a  variety  of  cases. 
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Section  7>  Tbs  tube  of  constant  croaa-aoetlon*  Performance 
curves ,  with  dlacuaalon^ 

In  Section  7,  which  In  many  waya  is  the  heart  of  this 

report,  we  diacuaa  the  curves  and  formulas  which  furnish  the 

ideal  performance  characteriatlca  of  the  hydropulse  for  the 

case  of  a  atraight  cylindrical  tube  where  any  of  the  parameters 

are  varied.   This  section  is  broken  into  five  parts  each 

devoted  to  a  special  quantity  or  to  a  special  aspect  of  our 

problem,  aa  follows:  a)  The  effective  velocity  I  /M  , 

b)  The  time  of  stroke  t„,  c)  The  average  thrust  P  =  I  /t  . 

"  s    a/  a' 

d)  Discussion  of  various  quantities,  such  aa  the  gaa  preaaure 

and  exit  velocity,  aa  functiona  of  the  time,  e)  Effect  of 

a  gradual  buildup  (in  time)  of  the  gaa  preaaure. 

a)  The  effective  velocity 

The  formula  for  the  effective  velocity  (cf.  equation 
(6.12)  of  the  preceding  aection)  la 


(7.1) 

with 
(7.2) 


/2R0 

TUf                ( 

.-.[^ 

fl       "^1 

_i_        1                   «VI#1 

r 

Pi  1.  ^ 

+  _  ,  and 

5r 

^e 


(7.3) 


For  ease   in  computation  of  the  effective  velocity  we 
give   three  grapha,  Figurea  7.1,   7.2,   7.3,  aa  follows: 
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Figure   7.1  la   a  graph  of  J  2 — =  |ji  as   a  function  of   the   "fuel 

parameter",   i.e.    the  ratio  M  /M^  of  the  mass  of  gas  per  stroke 
to   the  mass  of  water  in  the    tube   initially}     both  quantities 
are   plotted  on  a   logarithmic  scale,   for  a   series  of  values   of 
the  initial  temperature    0.   ranging  from  3000°  Kelvin  down 
to   500°  K.     Figure  7.2  gives   the  quantity  h  (cf.    (7.2))  as   a 
function  of  the  pressure  ratio  P«/t>Q  for  various  values  of  the 
adiabatic  exponent     7.     Figure  7«3  shows   the  quantity 

(Vh  +  k     -  k)   (of.    (7.1))  as  a  function  of  k  for  various  values 
of  h. 

Computation  of  the  effective  velocity  in  any  given  case 
is   then  carried  out  as  follows:      Compute  first  the  quantity 


i- 


^ 


1 


(dlmenslonally  a  velocity)  which  occurs  In  several 


'L....    ™..  _ U^ZK 


formulas.  The  value  of  1/2—=  «—  Is  then  fixed  by  the  graphs 

r^  g 

of  Figure  7.1  once    the  ratio  M  /M     is   known.      The   quantity  k 
Is   than  fixed  by  (7.3)  once   the   initial  velocity  U^    Is  given. 
For  a  given  pressure  ratio  P^/Pq   the  quantity  h  is   then 

determined  from  the  graphs  of  Figure  7.2,   after  which  the 

r~2,       2 
quantity  vh  +  k  -  k  can  be  determined  from  the   graphs  of 

Figure  7.3.     In  this  way  both  factors   in   the  formula    (7.1)  for 

the  effective  velocity  are  determined. 

Finally,  Figure  7.4  gives  graphs   showing  the  dependence 
of   the  effective  velocity  on   the   fuel  parameter  M    /M     for 
various   initial  velocities  IT,.      These   curves  are  all  drawn 
assuming  h  =  1.1,    /J  =  2,   and    ^^    =  2000°K. 

The  essential  information  regarding   the  dependence  of 
the  effective  velocity  on  all  parameters   is   clearly  contained 
In  these  four  sets  of  curves.      Certain  pertinent  general 
conclusions   should,   however,   be  drawn  from  them.      To  begin 
with,  we  observe   that  the  quantity  h  (see  Pigxrce  7.2)   Increases 
only  slightly  with  Increase  of  the  pressure  ratio  P^ /p     once 
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the  pressure  ratio  becomes   larger  than  30  or  40,      Consequently 

the  same    Is   true  of  the  effective  velocity,   since   it  depends 

on  the   pressure  ratio  only  through  the  quantity  h  (of.  equations 

(7.1)  and  (7.2 ))«     Prora  this  we  draw  an  important  practical 

conclusion:      The  efficiency  of  the  hydropulse  is  not  greatly 

increased  by  going  to  pressure  ratios  much  above  40.  We 

observe  also    that   the  values  of  h  for  any  reasonable  value 

of  7  and  for  pressure  ratios  near  40,   say,   are   close   to  unity. 

Figure  7.4,  which  is  drawn  for  h  =  1.1,   thus  furnishes  graphs 

for  reasonable  pressure  ratios.      Since   the   initial  velocity 

JJ.   of  the  water  occurs  only  in   the  quantity  k  in  (7.1),   it 

follows   that  Figure  7.3  indicates   the  behavior  of  the  effective 

velocity  with  increase   in  initial  velocity*      The  effective 

velocity  decreases  steadily  with  increase  in  TJ,,   tending  to 

zero  as  U,    increases  indefinitely.      This  was,  of  course,   to 

be  expected  since   the   thrust  augmentation  naturally  decreases 

as  U^  increases. 

Figure  7.4  is  one  of  the  moat  important  graphs,   since 

it  gives  an   indication  of  the  efficiency  as  a  function  of  the 

quantity  of  fuel  used  per  stroke.      The  effective  velocity 

decreases  steadily  with  increase   In  the  fuel  parameter,  with 

its  maximum  occurring  for  M  ^^     =  0.      The  maximum  Is   Infinite 

when  the  initial  velocity  TJ.   is  zero,  but  la  finite  for  non-zero 

initial  velocities  and  steadily  decreases  as  TJ^  increases. 

This  behavior  of  the  effective  velocity  comes  about,   of  course, 

because   the   thrust  augmenting  mass  M     grows   large  with  respect 

to  M_.      One  should  also  bear  in  mind   that  the  pressure  ratio 

S 
Pj,/P     is  considered  to  be  kept  constant  —  it  la  not  at  all 

obvious   that  that  could  be  accomplished  in  practice  when  M 

Is  very  small.     The  practical  conclusion  to  be  drawn  is   thlst 

The  quantity  of  fuel  injected  per  stroke  should  be  kept  as 

low  as  possible,   insofar  as   that  can  be  reconciled  with  the 

attainment  of  other  objectives   such  as  the  creation  of  a  pertain 

mlnlBium  average   thrust. 
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b)   The  time  of  stroke  tg 


The  formula  for  the  time  of  stroke  t_  (see  (6.13))  is 


(7.4)  tg 


(/)7  -   1 
^e 


s 


1.%  J^*  yt 


.2 


The   quantities  h  and  k  are  defined   in   (7.2)   and    (7.3). 

The  formula  for  t„  shows   clearly  that  t_   tends   to  zero 

S  3 

when  the  fuel  parameter  M   /^  tends   to  zero   (since  k-»-  oo  when 

jj    /M     -♦  0.  )     At  first  sight  this  seems  rather  paradoxical, 

since  one  would  be  inclined  to  argue  that   the   time  required 

for  an  entire   stroke  would   tend  to  grow  very  large  when  the 

amount  of  fuel  injected   is  very  small.     However,   it  must  be 

remembered  that  M     alone   is  being  varied  while   the   Initial 

pressure  p^    is  kept  fixed.     This  means   that  the  initial  volume 

V^  Is  being  made  smaller,   so  that  a  progressively  smaller 

volume  of  water  would  be  forced  out  of  the   tube  i)er  stroke, 

during  the  expansion  of  the  gas  down  to  the  exit  pressure  p  . 

Hence  one  sees   that  it  is  after  all  not  unreasonable   that 

tg  -••  0  when  M   -^  -#•  0  if  T?^/Vq   is  l»ld  fixed.      On  the  other 

hand,   when  M^L  grows   large,    it  is  readily  seen  that  t 

increases  essentially  as  the  square  root  of   this  quantity, 

since   then  k  tends  to  zero. 

As  an  aid  In  the   calculation  of  t„  from  (7.4)  we  have 

s 

drawn   (see  Piguire  7.5a  and  7.5b)   curves  for  the  third  factor 
on  the  right-hand  side  of  (7.4)  as  a  function  of  the  pressure 
ratio  P^/Pg  for  various  values  of  k;      two  such  sets  of  curves 
are  given  --  one  for     7  =  1.4  and  another  for     7  =  1.2. 
Prom  Figure  7.5  we  can  inf«r  how   t     varies  with  respect  to 
changes  in  the  pressure  ratio  and   the  initial  velocity.     As 
we   see,   t^  decreases  as   the  pressure  ratio  increases   (beyond 
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the  value  10,  at  least)  but  the  decrease  In  time  of  stroke 
Is  rather  slow  for  pressure  ratios  above  40.   The  time  of 
stroke  of  course  decreases  as  the  initial  velocity  U^  increases. 

c)   The  average  thrust  Pg  =  Ig/^a* 

^-~-^^— ^-"--^^^ ' — ■ ■ —  —  —    —         -     ^1 

The  fomula  for  the  average   thrust,   as  derived  In 
the  preceding  section   (of.    (6.14))  is 

(7.5)        Pg   =  :^  =  2p^Qg  ^1 Jh^+  k^    {Jl?+  k^-  k). 

^e 

The   quantities   h  and  k  are  defined   in   (7.2)   and    (7.3). 

Figures  7.6a,b  give   two  sets  of  graphs  of  P  /p^Q^  as 

a  function  of  P^/Pg  for  various  values  of  k   (i.e.   for  various 

Initial  velocities  U^ ).     For  one   set  of  graphs      7  =  1.4,   and 

for   the   other     7  =  1.2.      A   comparison  of„(7.5)  with  (7.3)   shows 

that  Fg  depends  upon  the  fuel  parameter  w^  and   the  initial 

w 
velocity  Uj   only  in  the   combination  given  in   the  formula   (7.3) 

for  the  quantity  k.     Figure  7.7  therefore   includes   the  effects 

of  both  fuel  parameter  and   initial  velocity  on  F     by  giving 

Fg/PgQg  as  a  function  of  k.      Two  sets  of  graphs  are  given J 

one  each  for     7  =  1.4  and      7  =  1.2  and  for   the  values 

P^/Pq   =  10,   40,    100. 

As  we  would   expect.  Figures   7.6a,   b  show   that  the   average 

thrust  P     increases  steadily  with  increase  of  the  pressure 
s 

ratio.      The  increase   is  rather  rapid  at  first.      But  for 

Pj^/Pq  >  40,   say,  Fg  increases  nearly  linearly  with  Pji/Pq* 

From  Figures  7.7a   and  7.7b  we   can  deduce   the  effect  of 

M 
changes   in  either  the  fuel  parameter  ^  or  the   initial  velocity 

w 
U,    on  F   ,   since   these   two  quantities  occur   in   (7.5)  only  throiigh 

the  quantity  k  given  by  (7.3).     As  ija  -►  0,   k  -*  00,  while 


w 
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k  -»•  0  when  -A.  -►  oo  •      The  quantity  k  is   linear  in  Uj^.      As 

we   see  from  pfgures  7.7e  and  7.7b,  Pg   steadily  decreases  from 

a  finite  value  at  k  =  0  and  approaches  asymptotically  another 

constant  non-zero  value    (which  can  be   readily  proved   from  7»5 

to  be   just  half  the  value  at  k  -  0)  as  k  tends   to  <X)  •     In 

other  words,    the  values  of  F.  lie  between  two  finite  bounds 
'  s 

(with  the  lower  bound  half  the  upper  bound)  for  all  values  of 

the   fuel  parameter  M  /M     and   the   initial  velocity  TJ,»      This 

is   a  rather  remarkable  result.      It  says,   for  example,    that 

the   quantity  of  fuel  injected  per  stroke  can  be  decreased 

indefinitely  and  the  average   thrust  over  the   stroke  will  not 

tend   to   zero,   but  will   tend  rather   to  a   limit  which  has   half 

the  value  of   the  maximum  which  is   approached  when  the   quantity 

of  fuel  injected   is  made   large  without  limit.        This  result 

is  made  somewhat  less   surprising  when  we  recall  that  both  the 

impulse   per  stroke  I_  and  the   time  of  stroke   t^   tend   to  zero 

when  M    /ll     tends   to   zero:      apparently   they  approach  zero   in 

such  a  way  that  the  ratio  I_/t„   tends  to  a  finite  and  non-zero 

limit  in  general.     It  should  also  be  recalled   that  the  pressure 

ratio  is  assumed  always   to  be  fixed  as   the   limit  is  approached, 

and,  as  in  a   similar  case  above,    it  is  far  from  clear  that 

M  /M^  can  be  made  very  small  while   Pj^/p^  is  maintained  fixed. 

In  other  words,    it  seems   likely   that   the   pressure  ratio  which 

can  be  achieved  on  Injection  of  fuel  will  tend   to  zero  with 

the   quantity  of  fuel  injected,   and  so   the   average   thrust  would 

also  tend  to  zero  through  failure   to  maintain  a  minimum 

pressure  ratio.     Nevertheless,   in  view  of   the  fact  that  th9 

efficiency  as  measured  by  the  effective  velocity  I    M     has  Its 

s/  g 

maximum  for  M  /M^  -  0,  the  following  conclusion  is  strongly 
indicated:   The  hydropulse  should  be  operated  at  high  frequency 


It  is  to  be  remembered  that  the  average  thrust  in  question 
here  refers  to  the  average  over  the  power  stroke  alone  and 
not  over  the  full  cycle. 
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ulth  Injection  of  spall  quantities  of  fuel  per  stroke*     Tha 
quantity  of  fuel  Injected   should  be  kept  as  loia  aa  la 
compatible  Yilth  maintenance  of  a  minimum  pressure  ratio   (of 
the  order  of  30  to  40  if  possible.) 

The  actual  numerical  values  taken  on  by  P-/P-Q-,   as 
indicated  in  Figures  7.7a,b  are   also  not  Twithout  interest. 
These  values  are  all  seen  to   lie  between  1»75  and  12,  depending 
on    y    and  the  initial  pressure  ratio   (which  we  have  assumed 
to  be  between   10  and   100.)     We  recall   (cf.    the  discussion  in 
Section  3  on   the  meaning  of   the   quantities  I  and  P)   that  the 
quantity  P     furnishes   the  average   thrust   (over  the  power  stroke 
only)  that  is  available  for  overcoming  the   total  drag  on   the 
hydropulse.      If  the  drag  is  given,   the  graphs  7.7a, b  make  it 
possible   to  estimate   the  minimum  cross-section  area  of  the 
hydropulse   that  would  be  necessary  to  overcome  it.     For  example, 
suppose   that  a  drag  of  1500  pounds  is  to  be  overcome  by  a 
hydropulse  which  is  operating  quite  near  the  surface  of  the 
water  so   that  p     =  15  pounds  per  square   inch.     Suppose,   in 
addition,   that  time  t     for  the  power  stroke  is  only  1/3  of 
the  full  cycle   (i.e.    the  acavenging  and  refill  parts  of  the 
cycle  require   twice  aa  much  time  aa   the  power  atroke.)     If 

7~  1.2  and  the  pressure  ratio  P^/Pg  is  about  40,    then 
Pg/PgQg  ~  ^  approx.  (cf.  Figure  7.7a).     The  average   thrust  Tg 
over  the  full  cycle  would  therefore  be  given  by  T„  A)^Q-   -  1. 
For  the   cross-section  area  Q     of  the  hydropulse^  100  square 
inches  would  then  be  requii^ed  for  a  thruat  of  1500  pounds. 
In  other  words,  a  cylindrical  tube  more   than  10  inches  in 
diameter  would  be  required. 

d)     Gas  pressure  and  exit  velocity  as  functions  of  the 
TEne. 

Our  approximate  foimulaa  for  the   impulse,   time  of 
stroke,  and  average   thrust,  as  diacuaaed  above,  were  all  derived 
on  the  basis  of  results  obtained   through  integration  of  the 
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equation  of  motion  of  ths  water  In  the  hydropulse.     For  this 
reason,   and  also  because   some  of  the  results  have  an  Interest 
for  their  own  sake,  we  discuss  briefly  the   integration  of  the 
equation  of  motion   (3.6)  for  the   case  of  the  tube  of  constant 
cross-section. 

In  terms  of  the  dimensionless  variables  introduced  in 
Section  4,   the  differential  equation  to  be   integrated  is 

(7.6)  v'^v-'y  -  1  -  (1-v)  ^  (v*)2  =  0  , 

since  in  this  case  the  quantity  q  =  Q/^g  is  identically  one. 
The  prime,  we  recall,  refers  to  differentiation  with  respect 
to  our  dimensionless  time  t.  Since  ^(v  )  =  2v  we  may 
write  (7.6)  in  the  form 

(7.7)  2(l-v)  v"  =  v'Vv"'^'  -1  . 

This  differential  equation  is  to  be  Integrated  for  the  initial 
conditions 

(7.8)  (     v(0)   =  ▼ 


j      v(0)   =  Y^ 


i.e.  for  a  given  initial  gas  volnnie  and  its  rate  of  ezpansioa 
The  latter  quantity  is  fixed  by  the  initial  velocity  of  the 
water  in  the  tube  at  the  moment  fuel  is  injected)     it  corre- 
sponds,  therefore,  roughly  to  the  forward  velocity  of  the 
hydropulse  since  the  water  in  the   tube  at  the  end  of  the  refill 
portion  of  the  cycle  would  have  about  this  velocity.     The 
quantity  Vm   la  fixed  by  the  initial  pressure  of  the  gas  once 
Its   temperature  and  molecular  weight  are  known.     Once  v  has 
been  determined  from  (7.7)  and  (7.8)  as  a  function  of   '''  ,  we 
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know   that  all  quantities  of  Interest,   such  as  the   Impulse, 
velocity  of  the  water,   etc.   can  be   computed  at  once* 

Our  problem  contains  foxor  parameters:     v     '^  V  /V  , 

S      S '   6 

that  Is,  the  ratio  of  the  gas  volume  at  pressure  p  to  the 
total  volume  of  the  tube,  7  the  adiabatic  exponent,  and 
V.  and  v^. 

The  Integrations  have  been  carried  out  numerically  In 
quite  a  large  number  of  cases  by  the  methods  discussed  in 
Appendix  II  of  this  report*  We  report  here  on  a  few  typical 
cases  in  the  form  of  graphs.  In  the  cases  we  report  the 
adiabatic  exponent  has  the  value  7  =  1.3.   The  quantity  Vj 
has  the  value  .0256  —  that  is,  the  initial  volume  of  the  gas 
is  chosen  somewhat  less  than  3  per  cent  of  the  total  volume 
of  the  tube.  The  initial  pressure  ratio  of  the  gas  Vt/V     is 
40  and  this,  together  with  the  adiabatic  gas  law  in  the  form 

(7.9)  pv-V  -  p^v];  =  p^v-V  , 

has  the  effect  of  fixing  the  quantity  v  in  (7.7). 

The  results  of  the  numerical  integration  are  shown  in 
the  graphs  of  Figures  7.8,  7.9,  and  7.10.  Figure  7.8  gives 
the  "excess  pressure  ratio"  V/P^-   1  as  a  function  of  the  "time" 
T  for  Vj  =  0,  1,  3. 

It  might  be  noted  that  the  thrust  P  at  any  time  is 
given  in  this  case  by  the  formula 

(7.10)  -  P  =  P^Qe^PAe"  ^^  ' 

ao  that  Figure  7,8  yields  P/i>gQQ  as  a  function  of      r.     The 

^  t 

If  p  is  atmospheric  pressure,  these  values  of  v^  correspond 
to  forward  velocities  of  the  hydropulse  of  0,  31.8,  and 
95.4  miles  per  hour. 
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important  observation  to  be  made  from  the  graph  is  that  the 
pressure  falls  off  with  extreme  rapidity  Initially  —  a  fact 
which  has  been  repeatedly  stated.  We  observe  also  that  the 
rate  of  fall  increases  with  increase  in  v^^,  that  is,  with 
increase  in  the  initial  velocity  of  the  water  in  the  tube. 
For  higher  pressure  ratio  also  the  rate  of  decline  increases 
(though  we  do  not  reproduce  here  our  results  for  higher 
pressure  ratios.)  Figure  7.8  Indicates  why  our  approximate 
formulas  based  on  the  assumption  of  Impulsive  action  of  the 
injected  gas  on  the  water  are  quite  accurate.  Figure  7.9 
makes  this  even  more  plausible.  In  these  graphs  the  rate  of 
expansion  of  the  gas  is  shown  as  a  function  of  the  gas  volume, 
which  is  the  same  effectively  as  giving  the  velocity  of  the 
water  as  a  fxinction  of  the  gas  volume.  As  we  see,  the  velocity 
goes  up  very  quickly  at  the  beginning  of  the  stroke  and  then 
remains  nearly  constant  over  the  remainder  of  the  stroke. 
The  rate  of  approach  to  a  nearly  constant  velocity  is  faster 
with  the  higher  initial  velocities.  From  these  observations 
on  the  behavior  of  the  gas  pressure  and  the  velocity  of  the 
water  in  the  early  part  of  the  stroke  we  see  that  the  accuracy 
of  our  approximate  formulas  is  better  the  higher  the  initial 
pressure  of  the  gas  and  the  initial  velocity  of  the  water  are. 

Figure  7.10  is  a  graph  showing  the  gas  volume  as  a 
function  of  the  time.  As  we  see,  the  gas  volume  is  nearly 
linear  in  the  time  except  near  r  =  0  where  the  curves  all 
have  high  cxirvatures.  The  method  of  numerical  integration 
of  the  equation  of  motion  explained  In  Appendix  II  is  baaed 
on  this  observed  behavior  of  v  as  a  function  of   r , 


In  Sections  9  and  10  of  this  report  two  possible  means 

of  smoothing  out  these  initial  sharp  variations  in  pressure 

are  discussed. 
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e)  Effect  of  linear  buildup  (Id  time)  of  the  gaa  pressure. 

In  all  but  this  part  of  tha  present  section  we 
assume  that  injection  of  the  fuel  leads  instantaneously  to  the 
initial  gaa  pressure.  We  shall  investigate  briefly  here  the 
motion  of  th©  water  in  th©  tube  which  results  when  the  gas 
pressure  is  assumed  to  increase  linearly  in  the  time  from  the 
value  p  to  the  value  p.«  In  this  way  we  attempt  to  simulate 
what  might  occur  in  practice  upon  injection  of  a  chemical 
into  the  tube.  We  obtain  graphs  from  which  the  values  of  v(r_) 
and  V  (t  )  can  be  obtained  once  the  pressures  p,  and  p   and 
the  "time"   7"  required  to  attain  the  pressure  p,  are  given. 
Clearly,  when  v(t  )  and  v  ( 7"  )  are  once  known  our  usual 
methods  and  results  for  the  case  of  instantaneous  creation 
of  the  initial  pressure  can  be  used» 

We  proceed  to  derive  the  appropriate  differential 
equation,  making  use  of  the  dlmensionless  variables  used  above 
(Section  7d).   In  view  of  (7,7)  and  (7,9)  we  have 


(7.11)  v"(l-v)  =|(P/Pe-  1)  , 

in  which  P/Pq   is  assumed   to  be  a  linear  function  of      r ,  as 
follows: 


(7.12)  p/p^  =  1  +  2cT    ,   0  <T<  T     , 


Thus   (7.11)   can  be  written 

(7.13)  r'u  -  v)  «  cT    . 

Aa  initial  conditions  we  have  t(0)  ^  0,   naturally,  and  we 
take  V  (0)  «=  0  also  (as  was  indicated  above.)  If  wa  were  to 
develop  the  solution  of  (7.13)  in  a  power  series  under  these 
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Initial  conditions  one  would  find  readily  that  the  series 

3 
begins  with  the  term  in   r  : 


(7.14)  ▼  =  a,T  +  ... 


.3 
'3' 


This  being  the  case  it  is  convenient  to  introduce  a  nen 
Independent  variable  y,  replacing   ''",  by  the  relation 

(7.15)  y  =^r^   . 

In  this  variable  the  differential  equation  (7.13)  is 

(7.16)  (3y  1-;  +  2  3l)  (1  -  v)  P  2  , 

dy      ' 

and  tbs  initial  conditions  become 

-  V  =  0  ,  and 


(7.17) 


dv  >  T 


at  y  =  0,   the  condition  on  ^  at  y  =  0  resulting  from  (7.16) 
at  y  -  0,   in  view  of   (7.14).     By  solving   this  one   initial  value 
problem  we  are  enabled  to  determine  v  and  ▼     as  functions  of 
y  and,   since   the  constant  c  does  not  occur  in  (7.16)  or   (7.17), 
we  also  obtain  at  one   stroke  the  solutions  of  (7.13)  for  any 
value  of  c. 

Suppose   that  the  pressure  ratio  P/P^  la  to  have   the 
value  P./Pq  at  a  specific  time       t^.     The  relation  (7.12) 
would  then  determine   the   constant  c: 

Pi 


-1 
(7.18)  c  = — ga- 


1^ 


o 


CONFIDENTIAL 


61 

CONFIDENTIAL 


Our  solution  of  (7.16)  nould  then  yield  through  (V.IS)  values 

for  V  and  v  at  '^  ~     "^q  ^^^   from  these,  in  turn,  our  usual 

methods  would  determine  the  motion  over  the  remainder  of  the 

3  troke • 

The  two  graphs  shown  in  Figure  7.11  were  drawn  to 

facilitate  such  calculations:   c  is  determined  from  (7.18) 

and  from  it  and   t     the  value  of  y.  is  fixed  by  (7.15). 

o  o     I 

Prom  the  graphs  the  values  of  v(t  )  and  v  ("L)  can  then  be 

read  off,* 

It  is  of  some   interest  to   consider  a  few  numerical 

cases.     We  assume      7=     1.4  and    t    =  0.1,      Two  different 

o 

pressure  ratios  P^/Pq  are   taken,   i.e.    P^/Pq  =  20  and  Pj^/Pq   =  40. 

That  is,    the  gas   pressure  ratio  builds  up  linearly  to  a  value 

of  20  or  of  40  in  the   "time"     r     =0.1.      The  following  table 

summarizes   the  results  obtained  for  the  entire   stroke,  which 

comprises  both  the   time  during  which  the   pressure  ratio  builds 

up  from  unity  to   its  maximum  and   the    time  during  which  the 

gas  expands   adiabatically  down   to   the   exit  pressure   p^.      Thus 

the  total  impulse    j„,    the   time  of  stroke    r  ,  and   the  dimension- 

s  s 

less  energy  c_  (imparted  to  the  water  by  the  expanding  gas) 
s 

are  made  up  of  two  contributions 

Ja  =  ^o  +  h   ' 

T      =    T      +    T 

a         0         1  » 


« 


It  should  be  recalled  that  these   curves  refer  only  to  the 
case   in  which  the  water  in   the   tube   is   initially  at  rest. 
There  would  be   no  difficulty  in  treating   the   cases  with 
U.    ^  0,  but  the   calculations  would  be   somewhat  lengthier 


than  those  for  U^^  =  0. 


CONFIDENTIAL 


62 


CONFIDENTIAL 


* 


distinguished  by  subscripts,  the  subscript  zero  referring 
to  the  period  of  linear  increase  in  pressure.   (The  quantity 


€      is  defined  by 


*.  = 


V^' 


where       E    is  given  in  (5.1).) 
s 


Linear  Pressure  Increase 


Pi 

Jo 

h 

Js=VJl 

■^o 

"^1 

8    s>     1 

20 
40 

.475 
.975 

.274 
.648 

.749 
1.623 

0.1 
0.1 

.158 
.256 

.258 
.356 

Pi 
Pe 

^0 

h 

3       O       X 

v(r) 

v^r^) 

^ 

20 
40 

.226 
.951 

.336 

1.696 

.562 
2.647 

.0158 
.0325 

.475 
.975 

.317 
1.300 

The  quantities  v(  r  )  and  v  (t-  )   are   the  diinensionless  gas 

volume  and  rate  of  expansion  at  the   time    ts     t  ,      The  quantity 

m     is   the  product  m  =  v-    •   P^/Pq  of  the  volume  v,   and  the 

maximum  pressure  ratio;      these  numbers  are   computed  in  order 

to  make  comparative  estimates  of  the  quantities  of  fuel 

necessary  in  these   cases  and  some  others   to  be  discussed  shortly. 

Clearly  the  quantities  m     are  proportional  to   the  masses  of 

S 
gas  present  if  p^  and  the   temperature   and  molecular  weight  of 


The  quantity  j^   is  given  by:    j^=  ^      f   (^  -Ddr     =  §  r,      • 

o       ® 
We  determine   j,   by  our  usual  methods  from  the   initial 

conditions  Vj=  v(t  ),Vj^  =  v  (t  ). 
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the  gases  are    the   same    in   the   cases   to  be    compared   (assuming, 

as   always.    Ideal  gases.)     We   observe    that   the    time     r     =  0,1 

assumed  by  us    turns  out   to   be  about  40  per   cent  of   the  entire 

stroke   In  one   case    (p^/Pq   ~  20)  and  about  30  per  cent  in  the 

other    (when  P^^/Pq   ~  40).      However,    the   Impulse  delivered  during 

the   time  of  increasing  pressure   is  a   considerably  larger 

fraction  of  the   total  impulse  over  the  fuel  stroke   --  it  is 

about  65  per   cent  when  Va/Vq   ~  20  and   about  60  per  cent  when 

p./p     =  40.     We   conclude!      If   the    time   required  for  the   gas 

to  reach  its  maximum  pressure  is  an  appreciable  fraction  of 

the  total  time  of  stroke,    this  period   is  of  importance  because 

the   impulse  delivered  then  is   likely  to  be  a   large  fraction 

of   the   total  impulse. 

We  proceed  to  make  a   comparison  with  a  few  cases   in 

which  the   initial  gas  pressure  is  assumed  to  reach  its  mazimum 

instantaneously.      On  what  basis   the   comparison  should  be  made 

is  not  entirely  clear.     We   have   chosen  to  make   a  comparison 

in  two  different  waysj      1)     We  assume   that  P^/Pq  and  v^^  =  v(t  ) 

have   the  values  given  in   the  above   table,  and  compute    J     and 

T-  by  our  usual  methods  for  instantaneous  rise  of  pressure, 
^  t 

with  the  additional  assximption  that  v,    =  0   (that  is,    the  water 

is   initially  at  rest)  at     "^  =  0.     If  we  assume   that  the 

temperature  and  molecular  weight  of  the  gases  are   the   same  at 

the  instant  when  their  pressures  are  maximum,   it  follows   that 

the  masses  of  gas  are   the   same   in   the   cases   to  be   compared. 


Pi/Pe 

^s 

h 

Linear 
Increase 

Instan- 
taneous 

Linear 
Increase 

Instan- 
taneous 

20 
40 

.258 
.356 

.410 
.318 

.749 
1.623 

.575 
1.281 

Pi/Pe 

«s 

"g 

Linear 
Increase 

Instan- 
taneous 

Linear 
Increase 

Instan- 
taneous 

20 
40 

.562 
2.647 

.336 
1.696 

.317 
1.300 

.317 
1.300 
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These  cslculatlons  shom  that  a  linear  buildup  of  pressure 
may  be  favorable  for  the  operation  of  the  hydropulse,  as 
compared  with  the  instantaneous  creation  of  the  peak  gas 
pressure.   The  impulse  delivered  per  stroke  for  a  given 
quantity  of  fuel  is  considerably  greater,  about  30  per  cent, 
when  the  pressxire  increases  gradually.  At  the  same  time  we 
note  the  marked  differences  in  the  time  of  stroke   r^. 

The  above  comparison  is  based  on  assuming  equality  of 
temperatures  at  the  instant  when  the  gas  pressures  are  at  a 
maximum.  But  when  the  gas  pressure  builds  up  gradually,  a 
considerable  amount  of  work  is  done  before  the  gas  attains 
its  peak  pressure,  and  consequently  the  assumption  of  equality 
of  gas  temperatures  may  be  unreasonable.  2)     In  the  second 
mode  of  comparison,  we  assume  the  same  values  for  P*/Pq  and 

s 

are  tabulated  below: 


for  the  total  energy   c„  imparted  to  the  water.  The  results 


Pl/Pe 

-^3 

h 

Linear 
Increase 

Instan- 
taneous 

Linear 
Increase 

Instan- 
taneous 

20 
40 

.258 
.356 

.261 
•  393 

.749 
1.623 

.740 
1.585 

Pi/Pe 

S 

^i 

Linear 
Increase 

Instan- 
taneous 

Linear 
Increase 

Instan- 
taneous 

20 
40 

.562 
2,647 

.562 
2.647 

.0158 
.0325 

.0264 
.0507 

It  appears  that  the  gradual  buildup  of  pressure  produces  almost 
the  same  total  Impulse  as  the  Instantaneous  case.  The 
approximate  formula 


=  V».^3 
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derived  in  Section  5  is  thus  found  to  be  accurate  even  in 
the  event  of  a  gradual  buildup  of  pressure.   (The  reason 
for  this  is  that  the  mass  of  water  in  tl»  tube  remains 
sensibly  constant  during  the  buildup.) 
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FIG.  7.2 
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FIG.  7.5a 
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FIG.  7.6a 
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FIG.  7.7b 


CONFIDENTIAL 


85 


S 


±a± 


A^cR/^(3E_Tl:teufi: 


m 


fe*i^ 


^teJ-IQNJ  6r 


< 


FIG.  7.8 
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FIG.  7.9 
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FIG.  7.10 
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FIG.  7.11 
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Section  8«      Tubes  with  eonloal  portions. 

We  have  seen  In  Section   5  that  the   impulse  per  stroke 
Is  not  sensitive   to  changes   In  shape  of  the  hydropulse   tube 
If   the   Initial  velocity  of  the   water  in   the    tube   is   not  great. 
Ifowever,   if  the  initial  velocity  of  the  water  is  high,   it  is 
possible   that  the   impulse  may  be   considerably  affected  by 
changes   in   shape  of   the    tube.      In   this   section  v/e   give   the 
results   of   calculations   for   the    impulse   and   time   of   stroke  for 
three  different  shapes   of  tubes  with  front  portions  made  up 
of  conical  pieces   joined   to  a  cylindrical  rear  portion.      These 
three   cases   are   compared  with  that  of  the   cylindrical  tube  of 
uniform  section.      The   calculations  in   some  of  the   cases  were 
carried  out  in  two  different  ways   —  once  by  numerical  inte- 
gration of  the  differential  equation  of  motion  and  once   through 
use  of   the  approximate   formulas  developed  in  Section   5,    so 
that  some   idea  of  the  accuracy  of   these  formulas  for   tubes 
of  non-uniform  section  is  obtained. 

We  make  use  of   the  dimensionless  variables  defined  in 
Section  4.      The  four   tubes  we   consider  are   then  defined  by 
giving   the  dimensionless   cross-section  area   q  as  a  function 
of   the  dimensionless  gas  volume  v,   as  follows   (the  accompanying 
sketches   show  the   tubes  in  the  correct  proportions): 

Case  a) 


a)       .  q  = 
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Case  t)) 


b) 


q  = 


(-  2v  +  2)2/^  ,  0  <  V  <  I 
1      »  i  <  V  <  1 


Case  c) 


c) 


q  = 


(2v  +  1)2/^  *  0  <  V  <^ 
(-  2v  +  2)2/3,1  <  V  <^ 

1     »  ^  <  ▼  <  1 


Caae  d)  (The  cylinder) 


d) 


q  s  1 
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We   have   calculated   the  values   of   j     and       r     in 

•'3  s 

each  of   these   cases  when        7  =  1.1,    P^/Pq   ~  74,   v.    =  .01. 

Two  different  values  for   the   initial  rate   of  expansion  of 

I  t  t 

the  gas   V.   were   taken:      v.    =  0  and  v-    =  3.      The   first  value 

corresponds   to  zero  initial  velocity  of   the  water  in  all  four 

cases.      The   second  value   corresponds   in   the   Case   d)   of   the 

straight   cylindrical  tube   to  a  velocity  of  95.4  mi.^/hr.    if 

the  exit  pressure   p     is   atmospheric  pressure.      In   the   other 

cases,   the  velocity  of  the  water  is  different  in  different 

parts  of   the   tube.     Also,    in  Cases  a)   and  b)   the  velocity  at 

the  entrance    initially  is  not   the   same   as   It  is    in  Cases   c) 

and  d),   but  is  greater  in  Case  a)   and    less   in  Case  b)   than 

in  Cases    c)   and  d)    since   the   cross-section  area   at  the  entrance 

is   less   in  Case   a)  and  greater  in  Case  b)   than  in  the  other 

two  cases.      Probably  a  more   realistic   comparison  for  Cases   a) 

and  b)  would  result  if  the   initial  expansion  rate  v.   were 

adjusted   so   that  the  entrance  velocities  would  be   the  same 

initially  in  all  cases   to  be   compared. 

The  following   table  summarizes   the  results  of  our 

calculations.      The   columns  headed  by  J_(approx.)  and       r 

s  s 

(approx.)  refer  to  values  obtained  by  use  of  the  approximate 
formulas  of  Section  5: 


The  approximate  value  of   r^  was  obtained  using  formula 

(5.19)  with  the  equality  sign  (converted  to  dimensionless 
form. ) 
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v;  =  o 

Case 

Js 

Jg(approx.) 

T 
S 

Tg(approx.) 

C3    a) 

P    ^^ 
c) 

1       1    d) 

1,361 
1.367 

1.366 

1.373 
1.373 

1.373 

.407 
.377 

.392 

.373 
.360 

.353 

v;  =  3 

Cn    a) 
d) 

.3013 
.3587 
.3303 
.3014 

.3053 
.3595 
.331 
.3019 

.139 
.159 
.154 
.150 

.138 
.158 
.154 
.148 

Upon  comparing  the  values  of  i  for  the  four  different 

shapes  of  tube  when  v.  =  0  (i.e.  when  the  initial  velocity  of 

the  water  is  zero)  we  see  that  they  are  practically  the  same 

In  all  cases.   This  is  in  accord  with  the  remarks  on  this 

point  in  Section  5,  where  it  was  stated  that  the  impulse  per 

stroke  is  insensitive  to  changes  in  shape  of  the  tube  at  low 

initial  velocities  of  the  water.  However,  when  v.  =  3,  some 

considerable  divergence  in  the  values  of  j  is  noted:   The 

impulse  is  20  per  cent  greater  for  the  tube  of  Case  b)  (with 

the  flared  out  front  section)  than  for  the  Case  d)  of  the 

straight  cylinder.   However,  as  we  have  already  remarked,  the 

1 
Cases  b)  and  d)  for  v.  =  3  may  not  be  strictly  comparable 

because  the  entrance  velocities  of  th^"  water  would  be  con- 
siderably different  initially.   This  criticism  does  not  apply 
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,  .  t 

to  tho   comparison  of  Cases   c)  and  d)  nhen  v.    =  3;     we  note 

that  j     is  10  per  cent  greater  for  Case  c)   than  for  the 
straight  tube  of  Case  d).      These   calculations   indicate   there- 
fore  that  the   impulse  per  stroke  may  be  appreciably  affected 
by  the   shape  of   the  tube  when  the   initial  velocity  of  the 
water   (which  is  roughly  the   same  as   the  forward  velocity  of 
the   hydropulse)   is   high.      It   is   indicated   that  a  flared-out 
front  section  may  have  advantages  from  the  point  of  view  of 
maximizing  the   impulse. 

We  see  that  the  values  of  both  j_  and        r   are  given 
quite  accurately  by  the  approximate  formulas  when  v-    -  3. 
The   accuracy  of   these  formulas   is  not   as  high  for     Vj,    =  0, 
especially  for        t  ,  but  the  error  is   still  not  unreasonably 
high.     We  observe  also   that  the  approximate   formulas  yield 
values  for  i     which  are   too  high  and  values  for        r    which 
are   too  low  --  thus  verifying  in   these   special  cases   the 
general  irule,   given  in  Section   5,    that  these   formulas,  yield 
upper  and  lower  bounds  for   j_  and     r    respectively. 

S  9 

It  is   clear  that  the  results  of  this   section  are 
concerned  with  too  restricted  a  range  of  parameter  values   to 
be  conclusive.     More  work  should  be  done   in   this  direction. 
It  should  also  be  pointed  out  that  the  assumptions  leading  to 
the  hydraulic  theory  used   throughout  this  report  are  probably 
too  crude   to  deal  accurately  with  tubes  of  rapidly  varying 
cross-section  area   —  the   flow  problems  become  more  acute   in 
these   cases. 


CONFIDENTIAL 


100 


CONFIDENTIAL 

Section  9.      Straight   tube  yiith  mass  and   spring   for  elaatlc 
energy  storage. 

As  we   have   seen   In  Section  7,    the  gas   pressure  decreases 
very  rapidly  from  Its   Initial  high  value   at  the  beginning  of 
the   stroke   while   the  exit  velocity  rises  very  rapidly.     It 
seems   likely   that  efficiency  would  be   promoted  by  smoothing 
out  these   initial  abrupt  variations  in  pressure  and  exit 
velocity.      In   this    section  we   investigate   the  use   of  a  mass 
and  spring  as  one  possible  method  of  achieving  a  smoother 
operation  of   the   hydropulse.      That  a  mass   and  a  mechanical 
spring  would  be  used   in   actual  practice   is   perhaps  doubtful  ~ 
our  results  simply  give  an  indication  of  the   possibilities 
inherent  in  provisions  for  elastic  storage   of  energy  during 
the   early  part  of   the   stroke.      In   the   section  following   this 
one,   another   less  flexible  but  perhaps  more   practical  method 
of  achieving    the   same  kind  of  results   is   considered. 

We  consider  the   case  of  a  cylindrical  tube  which  la 
fitted  at  the   end  by  a  piston  of  mass  M  backed  up  by  a  spring, 
as   indicated   in   the  following  figure: 


Spring 


K 


Pc 


Q- 
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The  fundamental  assumptions  regarding  the  behavior  of  the  water 
and  gas  to  the  right  of  the  piston  are  the  same  as  those  made 
In  Section  2.   Consequently  the  equation  of  motion  of  the  water 
in  the  tube  ia: 

(9.1)  -  1  WW  =i(p  -  p  ), 

in  which  W  is   the  volume  of  water  in  the   tube  at  any  instant. 
Equation    (9.1)  results  from  the   fundamental  differential 
equation    (3.6)  of  Section  3  when  it  is   assumed   that  the  tube 
has  a   constant  cross-section  area  Q     throughout.      The  equation 
of  motion  of   the   piston  ia  given   by 

(9.2)  -    M    (S   +  V)   =  -  ^  (v.-  W  -  V)   -  p   , 

In  which  V  is  the  volume  of  gas  at  any  instant,  V_  ia  the 
combined  volume  of  gas  and  water  when  the  spring  ia  in  its 
unstrained  position,  and  K  is  the  spring  constant.  In  addition 
to  (9.1)  and  (9.2)  we  used  the  relation  between  the  gas  pressure 
p  and  its  volume  V  which,  as  always,  ia  taken  as  the  adiabatio 

y 

relation  pV  =  const. 

We  introduce  the  dlmensionless  variables  w,  v,  and  r 
as  in  Section  4  and,  in  addition,  the  quantities 

(9-3)  a^=i- ,  «  =  ^  , 


in  place   of  the  mass  M  and  the  reciprocal  of  the   spring  constant 
K.      The  quantity  ip  is   the  ratio  of  the  mass  of  the  piston  to 
that  of  the  mass  of  water  required  to  fill  the  volume  V^. 
Equations   (9.1)  and   (9.2)  become   in  the  new  variables: 
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(9.4)  -  2wT(i"   s=  v^v"'^  -  1 

9 

(9.5)  -  2»//(v»"+  v")   =  i  (1  -  u   -  v)   -  v^v""^  , 


in  which  ths   primes  now  refer  to  differentiations  with  respect 

to   the   "dimensionless   time"     r    and  v_   is   the  value  of  v  when 

s 

the  pressure  p  of   the  gas  is  p^»     As   initial  conditions   (at 
the  "time"      r  =  0)  we   prescribe,   of  course,   the   initial  values 
of  V,  w  and   their  first  derivatives,   or,  what  is   the   same 
thing,   the   initial  position  and  velocity  of  the  piston  and 
the  gas-water  interface.      Our  initial  conditions   are,    then: 

V  (0)   =  Vj^  ,  w   (0)   =  Wj^   , 

(9.6) 

V  (0)   =  Vj^  ,   w    (0)   =  Wj^   , 

in  which  V*,  Vj,  w,,   and  w^;,   are  given  constants. 

Once   the   solutions  v(r)   and  w(t)   of    (9.4)   and    (9.5) 
satisfying   the   initial  conditions    (9.6)   have  been  obtained, 
the   thrust,   gas  pressure,   impulse,  exit  velocity  of   the  water, 
etc.    can  be   calculated.      Of  particular   interest  are    the   thrust 
P  and  impulse   per  stroke  Ig,  which  are  given   (in  dimensionless 
form)  tiy 

e  © 

(9.8)  J3  =  13/PVe//^  =-\    f   [  i(l  -  »  -  V)  H-  11  ar 
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The   first  formula   results   at  once   from   the   fact   that   the 
thrust  in   this   case  can  he   computed  from  the   strain  in   the 
spring.     The  upper  limit     r     on   the   integral  in   (9.8)  refers 
to   the    total  time   of   the   stroke,   which  is   taken  here,   as   in 
the   preceding   cases,   as   the    time  required  for  the  gas   to 
expand  down  to  the  exit  pressure  p  • 

The  integration  of  the  differential  equations  has  been 
carried   out  for  certain   special  initial  conditions  and   special 
values  for   the  parameters.        In  all  cases  we  have   taken 

7  =  1.5  (which  would  be  somewhat  too   high  in  practice,  but 
a  convenient  and  reasonable  value  for   the   purposes  we  have  in 
mind)   and  v     =  0.5,  which  means   that   the  gas  volume   at  pressure 
Pq   is   just  one-half  the  volume  of  the   tube  when   the   spring  is 
in   its  unstrained  position.     Except  for  one   case   to   be  discussed 
a   little   later,  we   take  always  as   initial  conditions  v.=  w,=  0, 
which  means    tiiat  the  piston  and   the  water  in  the   tube  are  at 
rest  initially,   and  v^  =  0.05,  w^^  =  0.95,   which  means    that  the 
initial  volume  of  gas   is   5  per  cent  of  the  volume   of  the  tube. 
With    7  =  1.5  one  readily  calculates   the  ratio  P^/P     of  the 
initial  gas   pressure   to  the  exit  pressure   from  the   adiabatlc 
relation  and   the  values  of  v_  and  v.;      the  result  is 

S  3. 

(9.9)  PiAe   ~  '^^•^^* 

For     y.  -«ve  have   taken  the  value        0.02,  and  for  ip  the  values 
0.1,   0.4,   0.5,   and   1.0. 


The  differential  equations   appear  not   to  be  integrable  by 
explicit  means.      The   results  given  here   were   obtained  by 
step-wise    integrations,  using   the  method   of  finite 
differences.      The  method  of   calculation     is  explained  in 
Appendix  I. 

This  means   that  the  exit  pressure  p     acting  on   the  piston 
would   cause   it   to    compress    the   gpring  an  amount  equal  to  2 
per  cent  of  the    length  of   the   tube,   as  one   can  readily  see 
(cf.    (9.3)). 
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Figure    (9.1)   illuatratea    in  a  general  way  what  occurs 
when  the   gas   is   "injected"   and   the  water  and  piston  are   set 
in  motion.      The  values  of  v,  w,   and  P/P-   as  functions  of   r 
are   shown  for   the    case       ip  =  0.1   (mass   of   piston:     mass  of 
water  filling   the  whole    tube   s=  0.1).      One   observes    that  w, 
representing   the  volume  of  water  in  the   tube,   steadily  decreases, 
as   it  should   since   the   water   is  always   subjected   to  a  resultant 
force   in   the  direction   of   the  exit.      The    curve   showing  the 
variation   in  v,   which  represents   the  gas   volume,    is   not  of  so 
simple  a   character.      One  observes   that  the   gas  volume  first 
increases  during  about  one-third  of  the   stroke  from  the  initial 
value  0.05  to  a  maximum,   then  decreases  during   the  second 
third  of   the   stroke    to   a  minimum,    after  which  it  rises    to  a 
value   slightly  greater  than  0.5  at  the  end  of   the  stroke. 
This  means   that  the  piston,  which  was   at  first  set  in  motion 
to   the    left  with  resulting   compression  of  the   spring,   was 
finally  decelerated  and  brought  to  rest,   after  which  it  reversed 
its  motion  and  gathered  enough  speed   to  cause   a   compression 
of  the  gas  in  spite  of  the  fact  that  the  gas-water  interface 
was  moving  always    toward   the  right.      This   behavior  of   the 
system  is  also  brought  out  by  the    curve   showing   the   pressure 
ratio  P/Pq  as  a  function  of   the   time:     At  first   the  gas  pressure 
drops  very  rapidly,   then  rises  again  when   the  gas   is    compressed 
by  the  action  of   the  piston,   and  finally  falls   to  the  exit 
value. 

Figures   (9.2)  and   (9.3)  give  the  pressure  ratio  and 
the  thrust    f    (cf.    (9.7))  as  functions  of    r   for  various  values 
of     ip  ,     For   the   larger  values       ip  =  0.5  and        i/^  =  1.0,   i.e. 
for  the   larger  piston  masses,    the   pressure  decreases   steadily, 
but  for   the   two   smaller  values   it  oscillates.      This  means 
simply  that  the  lighter  pistons  would  oscillate  rapidly  enough 
to   cause   a   compression  of   the  gas  during   the   time  of   the   stroke. 
It  is   of  some   Interest  to  observe    that   the   pressure   ratio   in 
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the   case        i/^  =  0»4   just  failed   to  attain   the  value  unity  near 

T  =  0.6.      This   suggests   that  for  some  value   of    ip  near  0.4 

there  vjould  be  a  discontinuous   change   in   the   value  for   the 

tl3ie  of  stroke       r    ,   when  it  is  regarded  as   a  function  of     ip  • 

s 

Figure  (9.3)  shows  the  variation  of  the  thrust  over  the  stroke. 
V/e  observe  that  the  thrust  actually  becomes  negative  over  part 
of  the  stroke  in  the  cases   ip  =   0.1  and   il^  =  0.4.   This 
means,  of  course,  that  the  "thrust"  on  the  hydro pulse  during 
this  part  of  the  stroke  would  actually  be  a  drag. 

Thus  we  see  that  the  time  variation  in  gas  pressure  and 
thrust  can  be  changed  quite  radically  by  varying  the  mass  of 
the  piston,  once  a  suitable  spring  has  been  chosen.   In 
particular,  we  observe  from  Figure  (9.3)  that  the  time  of  stroke 
seems  to  increase  with  decrease  in  mass  of  the  piston,  though 
this  effect  is  not  very  pronounced. 

The  graphs  shown  in  Figure  (9.4)  gives  the  thrust  as  a 
function  of  the  time  in  a  case  which  differs  considerably  from 
those  just  under  discussion.  In  this  case  we  suppose  that  the 
spring  is  compressed  initially,  and  that  the  piston  is  moving 
toward  the  right,  i.e.  the  compression  in  the  spring  is  being 
relaxed,  at  the  instant  when  the  gas  is  injected.  If  the 
initial  compression  in  the  spring,  the  velocity  of  the  piston, 
and  the  initial  pressure  of  the  injected  gas  are  properly 
chosen  one  might  expect  a  smoother  distribution  of  the  thrust 
over  the  stroke,  since  the  thrust  would  begin  with  a  positive 
value  (Instead  of  zero,  as  in  the  preceding  cases)  and  could  be 
maintained  positive  -throughout  the  stroke.   The  graphs  of 
Figure  (9.4)  bear  this  out  rather  strikingly.  These  graphs 


« 


The  value  k  =   0.02  which  characterizes  the  spring  in  our 
cases,  was  chosen  so  that  the  effects  we  seek  would  occur: 
It  is  clear  that  too  stiff  a  spring  would  furnish  essentially 
the  same  results  as  a  rigid  wall,  while  too  weak  a  spring 
could  not  supply  a  thrust  unless  very  large  displacements 
could  be  permitted. 
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resulted  from  the  following   choice  of  parameter  values; 

ip^  1«0  (a  rather  heavy  piston,   in  other  words),         xs  0.02, 
V     =  0.5,   V,    =  0,05,   w.    =  1.143    (v<+  w.    j^  1,   as  before,   since 

3  '1  1  11  I 

the   spring   is   initially  compressed),   v^  =  0,  vi^  =  -  0.3375  and 

-  0.75   (that  is,    two  different  initial  velocities  were   taken.) 

t 
One  observes   that   the   case   in  which  ^^    -  "  0.3375  yields   a 

particularly  smooth  th-rust  curve. 

Summing  up,  we  might  draw   the  following  general 

conclusions: 

1.  Provisions   for  elastic  energy  storage   can  result 
in   quite   considerable   changes   in   the  distribution  of    the  gas 
pressure   and    thrust  during   the    time   of   the   stroke. 

2.  The   time  of   stroke       t     may  be   a   rather  unstable 

s 

quantity.  As  we  have  seen  (in  the  case    i^=  0.4  above),  the 
time  of  stroke  may  change  discontinuously  at  certain  parameter 
values.   If  the  hydro pulse  were  to  be  provided  with  a  sensitive 
valve  mechanism  which  would  function  only  when  the  gas  pressure 
is  very  near  to  p  ,  it  could  happen  that  the  actual  operation 
of  the  hydro  pulse  would  be  erratic  if  it  were  to  take  place 
near  such  unstable  values  of  the  parameters. 

3.  If  gas  injection,  valve  operation,  and  the  functioning 
of  elastic  energy  storage  meclianisms  can  be  properly  timed, 

it  appears  possible  to  achieve  quite  smooth  distribution  of 
the  tlirust  (and  other  quantities  as  well)  over  the  stroke. 
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FIG .  9. 3 


CONFIDENTIAL 


111 


CONFIDENTIAL 


me 


115 


CONFIDENTIAL 

Section  10«  Spring  effect  obtained  from  a  layer  of  gas. 

In  the  preceding  section  we  have  seen  that  advantageous 
effects  on  the  operation  of  the  hydro pulse  might  be  achieved 
through  the  cushioning  effect  produced  by  use  of  a  mass  and  a 
mechanical  spring.  It  Is  not  inconceivable  that  a  mass  and 
spring  might  be  used  in  practice,  but  it  would  seem  more 
reasonable  to  try  to  achieve  similar  effects  through  use  of  a 
layer  (or  layers)  of  gas  as  a  "spring"  together  with  a  layer 
of  water  as  mass.   In  this  section  we  give  the  results  of 
calculations  for  the  case  of  a  cylindrical  tube  into  which 
two  layers  of  gas  are  "injected"  —  one  at  the  closed  end  of 
the  tube,  the  other  further  along  the  tube  toward  the  exit 
section  —  with  a  layer  of  water  between  the  two  layers  of  gas. 
The  situation  is  as  Indicated  in  the  figure: 


Closed 
una 


In  practice   it  is  possible,   even   probable,    that  many  such 
layers  might  co-exist  in   the   tube  at  any  instant.      However, 
it  is   probably  not  worthwhile    to  undertake   the  very 
laborious    calculations  which  would   be  required  for   such 
cases  until  more   specific  information  on   the  possible 
ranges   In  the  values  of  the  parameters   is  available. 
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We  make    the   same   fundamental  assiunptlons  as   in  Section  2 
regarding  the  behavior  of  the  gas  and  water  layers,   but  since 
we   have   to  deal  here  with  two  different  layers  of  both  gas 
and  water,   it  is  necessary  to  supplement  the  quantities  defined 
in  Section  2  with  a  few  others.     We   introduce   the  following 
volvime  ratios   (cf.    the  above  figure)  as   the  principle  dependent 
variables: 

(10.1)  ^e 

e 

in  which  V  is,  as  usual,  the  total  voliame  of  the  tube.  In 
e   ' 

addition,  we  introduce  the  quantity   a  =  W^ /V^  for  the  relative 
volume  of  the  layer  of  water  between  the  two  layers  of  gas. 
The  equations  of  motion  in  terms  of  our  dimensionless  variables 
are  than  obtained  for  each  of  the  layers  of  water  in  essentially 
the  same  manner  as  in  Section  3.  We  obtain  the  following 
differential  equations  for  the  quantities  ^    and  77  : 


(10.2)  A""^  -rilr]''^   -  2a^"  =  0 


s^      'a 


(10.3)    TT^r?"^  -  1  -  2(^"+  r7")(l  -  <f-77-a  )  =  0  , 

In  which  the  primes  mean  differentiations  with  respect  to  the 
dimensionless  time  r  and  the  quantities  ^„   and  t]     are  the 
relative  volumes  of  the  gas  layers  when  the  pressure  in  these 

layers  has  the  value  p  at  the  exit.  We  assume  the  gases  to 

7 
expand  in  accordance  with  the  adlabatic  relation  pV  =  const., 

with  -y  =  3/2. 

In  addition  to  the  differential  equations,  initial 
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conditions  for  the  position  and  velocity  of  the  masses  In  the 
system  must  be  iinposed.  In  our  numerical  calculations  we  have 
in  all  cases  assumed  the  two  masses  of  water  to  be  initially 
at  rest;  thus  we  have  as  initial  conditions 

(10.4)  ^'(0)  =  r)\o)     =  0  , 

(10.5)  ^(0)  =  .e^;  r?(0)  =77^, 

in  which   <f  .  and   77,  are  constants  giving  the  initial  values 
of  the  relative  volumes  of  the  two  gas  layers. 

Our  system  of  equations  contains  five  parameters: 
^■ft      "H  A  t       ^a>       Vat   and  a  •     We  do  not  propose  to  make 

X      J.      s      s 

an  exhaustive  study  involving  variations  of  all  of  these  parameters, 
since  our  purpose  here  is  simply  to  gain  a  qualitative  impression 
of  the  possibilities.   Thus  we  limit  ourselves  in  what  follows 
to  the  cases  in  which  ^ „   =       ^„  =  «  =  x  *  that  is,  we  assume 

S  S  ^  T 

that  the  volumes  of   the  gas   layers  at  the  pressure  p^  are  ^ 
of   tte  volvune  of  the   tube,   and   that  the   layer  of  water  between 
the   two   gas   layers  has   the   same  volxime.      However,   various 
combinations   of  values   of   the   initial  relative   volumes       <^  ^^ 
and      77.    have  been   taken,  which  amounts  to  saying  that  different 
combinations  of   the   initial  pressure   ratios   P^/Pq   in  each  gas 
layer  were   taken.      The  seven  different  combinations  of  values 
used  in   the  numerical  calculations  are   indicated  in  the 
following   table: 

(f.    =  .025,    .01     with  no  second  layer 

<f,    =  .025,     r]^   =  .025,    .01 

(f.    =  .01   ,     r)      =  .025,    .0156,    .01. 

The  relations  between  these  values  for  the  volume  ratios  and 
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the  corresponding  values  of  the  pressure  ratios  P/P-  is  as 
follows: 

Voltune  ratios    Pressure  ratios 

•025  32 

.0156  64 

.01  125 

The  results  obtained  through  numerical  integration  of 
the  differential  equations  for  all  of  these  cases  are  given 
in  the  form  of  graphs.  Figures  (10.1),  (10.2),  (10.5),  and 
(10.4)  give  the  variation  in  pressure  at  the  closed  end  of  the 
tube  as  a  function  of  time  for  each  of  the  seven  different 
combinations  of  initial  pressure  ratios.  The  pressure  ratios 
are  plotted  on  a  logarithmic  scale,  while  the  "time"  r  is  also 
plotted  logarithnically,  on  two  of  the  sets  of  curves.  The 
principal  observation  to  be  made  from  these  curves  is  that 
considerable  modifications  in  the  shape  of  the  early  part  of 
the  pressure -time  curves  can  be  achieved  by  taking  various 
combinations  of  initial  pressure  ratios  in  the  two  layers  of 
gas.   Consider,  for  example,  the  combination  of  pressure  ratios 
32  -  125.  In  this  case  we  observe  that  the  pressure  at  the 
closed  end  of  the  tube  rises  initially  (instead  of  falling  very 
rapidly,  as  it  does  in  most  cases)  because  of  the  high  presstire 
in  the  second  gas  layer.  However,  the  peak  pressure  ratio  is 
about  92,  that  is,  considerably  less  than  the  initial  high 
value  in  the  second  gas  layer.  This  shows  clearly  that  the 
smoothing  out  of  the  thrust  curve  in  the  early  part  of  the 
stroke  is  feasible  through  use  of  a  cushioning  gas  layer.   The 
case  125  -  32  is  also  interesting:  here  the  pressure  curve 
drops  very  rapidly  at  first,  but  reaches  a  minimum,  goes  up 

The  method  of  finite  differences  was  used.  For  a  full 
discussion  see  Appendix  L  of  this  report. 
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to  a  maximum  and  then  falls  steadily.  This  means  that  the 
second  layer  of  gas  iias  at  too  low  a  pressure  Initially  to 
prevent  the  first  layer  from  expanding  rapidly.  Instead,  the 
second  layer  was  compressed  at  first  but  later  gave  up  some 
of  its  stored  energy  to  compress  the  first  layer  slightly. 
One  sees  from  thsse  examples  that  the  interactions  between 
the  two  layers  can  be  rather  complicated. 

Figure  (10.5)  is  a  graph  of  the  (dimensionless )  total 
impulse  delivered  to  the  tube  up  to  any  time  r  ,  as  a  function 

T 

of   T.  That  is,  it  is  a  graph  of  J  =  ^  J  (p/p^-  Ddr 

o 
against  r  ,  with  P/Pq  representing  the  pressure  ratio  of  the 

gas  layer  at  the  closed  end  of  the  tube.   One  observes  that 

the  greater  part  of  the  impulse  is  delivered  during  the  early 

part  of  the  stroke.   However,  comparison  of  the  case  32-0 

with  32  -  32  and  of  the  case  125  -  0  with  125  -  32  indicates 

that  the  introduction  of  the  second  layer  has  a  slight  smoothing 

effect. 

It  is  found,  as  might  have  been  expected,  that  the  time 

of  stroke   r  (that  is,  the  time  required  for  tha   gas  at  the 

s 

closed  end  to  expand  down  to  exit  pressure)  is  increased  by 
the  introduction  of  a  second  gas  layer,  the  increase  being 
smaller  as  the  pressure  of  the  second  layer  goes  up.   The  total 
impulse  delivered  over  the  stroke  is  also  increased  by  the 
introduction  of  a  second  gas  layer.   However,  the  average 
thrust  j/Tg  over  the  stroke  does  not  always  increase  when  a 
second  gas  layer  is  present.   Our  calculations  indicate,  in 
fact,  that  the  average  thrust  may  be  decreased  by  the  presence 
of  a  cushioning  layer  if  its  initial  pressure  is  sufficiently 
low.  These  results  are  summarized  in  the  following  table: 
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Initial 

Pressure 

Ratios 

^s 

J 

Vts 

1 

™Puel 
(based 

Effective  Vel- 
ocity Ratio 

on  32-0  as  unity) 

125  -       0 

0.205 

1.34 

6.54 

0.91 

125  -     32 

0.257 

1.59 

6.19 

0.66 

125  -     64 

0.273 

1.69 

6.19 

0.64 

125  -  125 

0.254 

1.77 

6.97 

0.60 

32   -        0 

0.290 

0.93 

3.21 

1.00 

32   -     32 

0.370 

1.29 

3.48 

0.70 

32   -  125 

0.342 

1.57 

4.59 

0.65 

The  last  colxann  was  obtained  by  assuiulng  the  gases  in  both 
layers   to  have   the   same  Initial  temperature  and  molecular 
weight  in  all  cases.      These  ratios  are  baaed  on  assuming   the 
case  32  -  0  to  be   the   standard,    the  other  numbers  representing 
the  relative  effective  velocities   (for  both  layers   combined) 
in  the  other  cases.      The  effective  velocity  increases  as   the 
mass  of  gas  increases   (as  in  Section  7.) 

Figures   (10.6)  and   (10.7)  are   the  results  of  an  attempt 
to  compare   the  effects  of  a  mass  and  a  mechanical  spring  with 
those  for  a  layer  of  gas  used  as  a  "spring".     On  what  basis 
one  should  make   the   comparison  is  neither  uniquely  defined  nor 
entirely  clear.     We  have  chosen  to  make   the   comparison  assuming 
the  masses  of  the  piston  and   the  water  layer  to  be   the  same 
(that  is,       ip  =  0.25  for  the   spring  case),  while   the   spring 
constant  was   chosen   to   correspond   to   the   "stiffness"   of  the 
gas  layer  at  a  pressure  ratio  of  31.6.     In  Figure   (10.6)   the 
pressure  on  the  closed  end  of   the   tube   is  given  as  a  function 
of    T  for  various  initial  conditions,  with  pressure   in  the 
case  of   the  spring  referring   to  spring  thinist  per  unit  cross- 
section  area.      The  nxunbers  ,  (125  -  125)  and   (31.6  -  125) 
attached   to  the  curves  referring   to   the  mechanical  spring  cases 
mean  that  the   spring  was   initially  compressed  an  amount  which 
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would  correspond  to  Initial  pressure  ratios  of  125  and  31«6 
at  the   closed  end.      One   observes    that  tte   curves  for   the 
comparable   cases   coincide  fairly   closely  in   the   initial  stages 
of  tbs  motion,   but  differ  widely  during   the   last  70  -  80  per 
cent  of  the   stroke.      The  curves  for  the  mechanical  spring  cases 
show  rather  heavy  oscillations  in  pressure    (the   thrust  becomes 
strongly  negative),  while    those   for   the   two  gas   layers  fall 
off   slowly  and  steadily   to   the  value  unity.      This  brings  out 
an  essential  distinction  between   the   two  types  of  cushioning 
devices:      the  mechanical  spring  Is  able   to  exert  a  backward 
pull  on  the   end  of  the  duct,  even  though  the  gas  pressure 
ratio   is  well  above  unity,   while   such  an  action   is   impossible 
with  a  gas  used  as   spring.      Also,   a  gas   layer  would   be   likely 
to  be  always  a  rather   "stiff"   spring  —  for  the  present  purix)se 
this   "stiffness"  was  assximed   to  be,  on   the  average,   about  four 
times  as  large  as  in  the   cases   treated  in   the  preceding  section 
(in  which     «  had   the  value   0.02.) 

In  Figure    (10.7)   the  pressure  in  the   second  layer  of 
gas   is  compared  with  the  gas  pressure  on  the  piston  in  the 
mechanical  spring  case.     As  one   sees,    the   curves   in  the   cases 
to  be   compared  follow  each  other  very  closely.     In  the  cases 
with  pressure  ratios  of  31.6  and   125  even  the   slight 
Irregularities  near        t  ^  o.l  occur  for  both  mechanical  spring 
and  gas  layer. 
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FIG.  10.5 
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FIG  .10.6  -    CONFIDENTIAL 
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FIG  .  10.7  CONFIDENTIAL 
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Section  11.      The  valve  less  hydropulae."^ 

Except  in  this  and  the  following  section  we  assume 
always   that   the   left  end  of  the  hydropulse    tube   is  closed  by 
valves  of  some  kind.     In  this  section  we   consider  the  possibility 
of  operating  a  hydropulse  without  valves. 

We   shall  study   the   simple   device   indicated   in   the   figure: 


We   assume   that  a   layer  of  gas  is   injected   into   the   tube  while 
water  is  flowing   through  it,    toward  the  right,   as  Indicated 
in   the  figure.     A  constriction  of   the    tube   is   provided  in  front 
of  the  gas  layer,   so   that  a   thrust  becomes   possible.      Our 
object  is  to   study  the  motion  during   the   time  required  for 
the  gas   to  expand  down   to   the  outside  pressure  p  .     In  other 
words,   we  consider  here  only  what  occurs  during  a  single   stroke. 


•K-  The   idea  of  treating   this  problem  was  suggested  to  us  by 
Ensign  L.    B.   Edelman,    of   the  Bureau  of  Aeronautics   Project, 
Annapolis,   Md. 

^'"*The   hydropulse  without  valves   should  be  distinguished  from 
the   hydroduct    (P.   Zwicky's   terminology.)      in   the  hydroduct 
the  gas  would  be    introduced  continuously,   not   intermittently, 
so   that  a   steady  flow  would  result. 
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We  make  the  same  assumptions  regarding  the  behavior  of 

the  water  and  gas  layers  as  in  Section  2.  We  denote  the 

volume  of  the  wide  section  of  the  tube,  i.e.  that  extending 

from  Xn  to  the  exit  section,  by  V.  and  introduce  the  following 
X  e 

notation  for  the  relative  volumes  of  the  various  water  and 
gas  layers: 

„  -  ^^-^5^^e 

(11.1)         «  o  e 

(X3-X2)Q^ 
V=  ^^. 

We  denote  by  q  the  ratio  of  the   cross-section  areas  at  the 
front  and  exit  sections: 


(11.2)  q  ==  Q^/Q^ 


• 


The  equations  of  motion  are  then  readily  derived, 
essentially  as  in  Section  3,  with  the  result: 

(11.3)  (1   -   Wg-   Vlg)"^   V"^-   1   =   -   2W3   Wg   , 

(11.4)        (1  -  Wg-  ^^r^vl'   1  +  i^2^^{l   -  1/q^) 

+  2W2(W2+  Wj^/q^)  =  0. 

The  dependent  variables  are  «„  and  v:  v^  ±3   a  constant  and 

w,  ia  determined  at  once  when  Wg  and  v  are  known,  since 

W2+  v  +  WjS  1.   The  quantity  v  is  the  value  of  v  when  the 

pressure  in  the  gas  layer  has  the  value  p^,  and  7  is  the 

-V     -V  ® 
adiabatic  exponent;  that  is,  p_v'  =  pv'  .   The  primes  refer 

to  differentiations  with  respect  to  the  dimenslonless  time  t' 

The  initial  conditions  are: 
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(11.5)  ^2'r=o  ~  "21   *  ^'rs=0  ""  ^j 


I  •  t  I 


^2^r^  -  ^21   '^    It=0  ~  ^1 


t 


that  is,  we  prescribe   In  effect  the   Initial  positions  and 
velocities   of   the  gas-"aater  interfaces. 

The  number  of  parameters   involved  here   is  rather  large; 
we   have  been  forced   to  restrict  ourselves    to  a  fev9  of   the  many 
possibilities.      In  all  cases  we   have   taken  v_   =  0.80, 
Wg.    =  0.20,   Vj    =  0,   q  =  0.25,   and       7  =  1.5.      This  means   that 
the  expanded  volume  of   the  gas   (when  at  pressure  p   )   and  the 
initial  volume   of   the  water  layer  extending   from  the    constriction 
to  the  gas   layer  are  each  1/5  of   the  volume  of  the  wider 
portion  of  the   tube.      The  assumption  that  v.    =  0  means   that  the 
initial  velocities  of  the   two  gas-water  interfaces  are  assumed 
to  be   the  same.     In  addition  the   cross-section  areas   of   the 
two  different  portions  of  the  tube  are  assumed  to  be   in  the 
ratio  114.     With  respect  to  the  remaining  parameters  we  have 
considered  the  following  four  cases: 

1)  Vj,  =  .008   (or,   P^/Pq   =  125),   n^^  =       0        ,  Wj^  =  0.20 

2)  Vj^  =  .008   (or,   Pj^/Pg   =  125),  Wg^^  =  0.279   »  »]l  ~  ^'^^ 

3)  v^  =  .008        (Pi/Pe    =  125)      ,   Wg^  =  0.279   ,  w^^   =  0.025 

4)  Vj^  =  .020     (P^/Pq   =  31.6)      ,   Wg^  =  0.279   ,  w^  =  0.025 

The   initial  pressure  ratios  P^/Pq  result  from  the  adiabatlc 


*  The  dimensionleas  velocity  of  the  left  interface,   namely 
A-  (■   y    ),    Is  equal  to  WgJ      that  of   the  right  interface 


I 
is   -« 


'3- 
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relation  (with  7  =  1»5).  In  Cases  1)  and  2)  the  volume  of 
the  constricted  forward  part  of  the  tube  la  1/5  that  of  the 
wider  part,  and  1/40  in  Cases  3)  and  4).  The  value  Wg.  =  0.279 
means  that  water  would  be  flowing  into  the  entrance  section 
of  the  tube  at  a  rate  of  about  40  ml»/hr.  if  p  were  1  atm. 
In  Case  1)  the  initial  velocity  is  zero.  Two  different  initial 
pressure  ratios  are  also  taken. 

In  each  of   the  above   cases,   the   solutions  of   (11.3)  and 
(11.4)  were  obtained  by  step-wise  numerical  integration   —  a 
simple,  but  laborious,   process.      The  results  are  given  in   the 
form  of  graphs  in  Figures   (11.1),    (11.2),   and   (11.3).     Figure 
(11.1)   shows   the  variation  in  the    (dlmensionless)   Impulse   j, 
given  by 

rT 


1    f    (i^  - 


(11.6)  J*'^^"!     1      (^  -  1)  dT  ,    (1  -  q) 

o         * 

The  quantity  p,    la    the  pressure   at  the   shoulder.      The   curves 
for  each  of  our  four  cases  are  shown,   together  with  two  cases 
of   the   hydro pulse  with  a  valve    (i.e.   with  the   tube   closed  off) 
located  at  the  position  where   the   left  hand  gas-water  interface 
is   located  initially  in  the   case   of  the  valve less  hydropulse. 
The   curves  for  Cases  1)  and  2)  fall  too  close   together  to  be 
distinguishable  on  our  graph,  which  means   that  the   impulse  ia 
hardly  affected  at  all  by  a  quite    considerable   change   in  initial 
velocity  of  the  water.      The   corresponding   curve  for   the  hydro- 
pulse  with  a  valve    (and   initial  velocity  of   the  water   the   same 
as   in  Case   2))   shows   the    total  impulse   to  be   greater  with  a 
valve.      A  considerable   decrease    in   length  of   the  narrower  forward 
part  of  the   tube    (in   the  ratio  of  8   :    1)   leads   to  a   sharp  loss 
of   Impulse,   as  we   see  by  comparing   the  curve   for  Case   3)  with 
that  for   the   Cases   1)   and  2).      This   is,   of   course,   what  one 

See  Appendix  I  for   the  details. 
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expects:      some   Impulae   is   lost  because   of   the    transport  of 
momentum  out  of  the  front  end  of   the   tube.      The   ctirve  for  Case  4), 
in  which  the    initial  pressure  ratio  is  31,6    (instead   of  125, 
as   in  the  other  cases)  and  in  which  also  the  narrow  front 
section  is  short,   indicates  a  considerable   loss  in  impulae  by 
comparison  with  the  curve  for  the   corresponding  case  of  a  hydro- 
pulse   closed  off  at  the   left  end. 

Pigtire    (11.2)  gives  graphs   of  the   pressure  ratio   in   the 
gas   layer  and  at   the   shoulder  where   the   tube  broadens   out. 
The   latter  quantity  Pi/Pg  was  determined  from  the  relation 

(11.7)  P3^/p^   =  p/p^    +  2W2  Wg. 

The  pressure  on  the  shoulder   (which  of  course   produces   the 
thrust)   thus  will  be  nearly  that  of  the  gas   if  Wg  is   small  — 
that  is,  if  the  acceleration  of  the  water  in  front  of   the 
shoulder  is   small.     The  graphs  of  Figure   (11.2)   indicate   that 
the   two  pressures  do  not  differ  greatly,  except  at  the  beginning 
of   the   stroke   in  Case   3),   in  which  the  "nose"   is  short  and   the 
initial  gas   pressure   high. 

It  is,    in  fact,    possible    to  draw  a  few  generally  valid 
(and  physically  quite   reasonable)   conclusions  from   (11.7)  and 
(11.4),  with  regard  to  the  effect  .of  varying   the  front  or  "nose" 
of  the   tube.     We   see   that  the   coefficient  of  Wg  in   (11.4) 
becomes   large  when  w,    is  made   large  or  q  is  made   small,   that 
is,  when  the  relative  volume   of   the  nose   is  made   large  or  its 
relative   cross-section  area  is  made   small.     Both  of   these 
changes,    in  other  words,   lead   to   large   values  of   the   "effective 
inertia"  of   the  water  in   the  front  part  of  the   tube.     In   these 
cases,   then,  Wg  can  be  expected   to  be   small.     Thus  a   long, 
narrow  nose  will  in  general  lead   to  a  high  thrust,   since  then 
p,/p     will  differ  only  slightly  from  p/p^   and,   in  addition, 
the  shoulder  area  on  which  p,   acts   to  create   the   thrust  is 
larger. 

COl^IDENTIAL 


142 


CONFIDENTIAL 

In  Figure    (11.3)   the  velocities  of   the  gaa-uater 
interfaces  are  shown.      The  results  are  what  one  expects.      The 
sharper  variations   in  velocity  are   correlated  with  the  higher 
Initial  gas   pressures  and  shorter  length  of  nose.     We  observe 
also  that  injection  of   the  gas  at  a  pressure  ratio  of  125  is 
just  sufficient  to  bring  the  water  in   the  front  part  of  the 
tube    (which  has   initially  a  velocity  of   the  order  of  40  mi./hr. ) 
to  rest. 

Perhaps   the  main  conclusion  to  be  drawn  from  the   analysis 
of  this   section  is   the  following:     As  far  as  the   total  impulse 
for  a   single   stroke   is   concerned,   we    see    that  a  hydropulse 
without  valves   could  be  designed   so   that  its  performance  would 
compare   quite  favorably  with  that  of   the  valved  hydropulse. 
To  do  this,   however,   a  fairly  long  narrow  forward  or  nose,  section 
would  seem  to  be  needed.     As  we   shall  see   in  the  next  section, 
such  a  design  is  likely  to  cause  grave  difficulties   in  the 
"scavenging"   part  of   the  cycle,  during  which  the   tube  should 
be  refilled  after  the  power  stroke. 
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Section  12«   Scavenging* 

In  this  report  we  have  in  general  concentrated  our 
attention  on  a  single  exhaust  stroke  of  the  hydropulse  and  have 
left  aside  problems  concerning  other  phases  of  the  full  cycle. 
In  particular,  the  phase  of  "scavenging."  By  this  Is  meant 
the  portion  of  the  cycle  during  which  the  tube  Is  being  refilled 
with  water  In  preparation  for  the  Injection  of  fuel  and  the 
beginning  of  a  new  exhaust  stroke.  In  the  present  section  our 
aim  la  to  make  a  study  of  the  scavenging  process  In  .the  simplest 
possible  case  In  order  to  get  at  least  a  qualitative  Idea  of 
what  to  expect.  As  we  shall  see,  the  results  Indicate  that 
the  scavenging  phase  of  the  cycle  may  well  present  problems  of 
crucial  importance  for  the  hydropulse  development  as  a  whole. 

We  consider  the  simplest  possible  model  of  a  hydropulse 
with  a  valve:  a  straight  tube  with  a  bank  of  valves  inside 
the  tube,  as  indicated  in  the  figure. 


Q 


Intake 


L.Khau5t    Charr.ber 


Exhaust     Flow 
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The   problem  we  wish  to   consider  is    the  following:     We   assume 

that  the  exhaust  stroke  has   just  been  completed,   i.e.    that  the 

gas  has  expanded   (with  the  valves   closed)  until  its  pressure 

has  decreased   to   the  value  p_   on   the  outside  of  the   tube  and 

the  water  in   the  exhaust  chamber  has  acquired  a   certain  velocity. 

An  instant  later  the  gas  will  be  at  a  pressure   less   than  p  , 

(since    the   inertia  of   the  moving  mass   of  water  will   lead   to  a 

further  expansion  of  the   gas),    so   that  valves   could  be   opened 

automatically  by  the   pressure  of   the  water  in  the  forward  part 

of   the   tube.      Our  purpose   in  this   section  is   to  atudy^the 

subsequent  motion  of  both  layers   of  water,    on   the  basis   of  our 

usual  assiamptions  regarding   the   behavior  of  gas   and  water  layers. 

In  particular,  we   assume   the  exterior  pressure   p^  to  remain 

at  the   same  constant  value  at  both  ends  of  the   tube   throughout 

the  motion.     In  addition,  we  assume   that  the  valves  have  a 

negligible  effect  on  the  flow  once   they  have  been  opened. 

We  denote  by  V     the    total  volume   of   the   tube   in  the 

"exhaust  chamber"   to  the  right  of   the  valves.      The  quantities 

W,   and  Wo  represent  the  volumes  of  water  at  any  instant  in  the 

"intake"  and  "exhaust"   sections  of  the   tube  respectively   (see 

the  figure),  while  W,  represents   the  volume  of  gas  at  any 

instant.      These  quantities  are   therefore  all  functions  of  the 

time.     It  is  also  convenient  to  introduce   the  quantity  W     for 

the  volume  of   the   tube   to  the  left  (I.e.   in  front  of)   the 

valves,  which  means   that  W^  is   the  value  of  W,    at  the   initial 
'  o  1 

instant.     The  equations  of  motion  are  obtained,   in  the  usual 
way,    in   terms   of   the  relative  volumes  w,    *=  W./V   ,  with  i 
ranging  from  0  to   three.     They  are: 

n 


(12.1) 


S'l^l  =  1   -   (1  +  %-  'i-  '2^"^  '3a 


n 


SWgWg    r=    1    -    (1    +   w^-    »3_-    Wg)-^    Wg^g 
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in  which  the   primes  mean,   as   always,   differentiations  with 
respect  to  the  dimensionless   time     t  and  the  relative  voltane  w, 
of  the  gas  layer  has  been  replaced  by 


(12.2) 


w,   =  1  +  w^- 
o  o 


'l-''2 


w,-   is    the  value   of  w,  when   p  -  p  •      Since  w     is  a   constant, 
the   two  equations    (12.1)   serve    to  determine   the    two  quantities 
w,   and  Wg,  once  appropriate   initial  conditions  have  been 
prescribeci. 

We  report  here   the  results  of  integrating   (12.1)   in  the 
following   special   case:     We   assxAme   the   volume  of  the    tube   in 
front  of  the  valves  to  be   the  same  as  w^-*    that  of  the  gas 
layer  when  the   latter  is  at  the  exterior  pressure  p  ,  while 
the  exhaust  chamber  volume  V.   is  four  times  as  great.      (In 
other  words,    the  exhaust  chamber  is  four  times  as  long  as   the 
intake   section.)     Initially,    the  water  in  front  of  the  valves 
is  at  rest,  while   that  in  the  exhaust  chamber  has   the  velocity 
acquired  during  the  previous  expansion  of  the  gas  layer.      The 
latter  initial  velocity  was  obtained  from  previous  calculations 
(discussed  at  the  end  of  Section  7)  on  the  assumption  of  an 
initial  pressure  ratio  P-/Pq  in  the  gas  layer  of  125  and  1.5 
as   the  value  for  the  adiabatic  exponent     7  •     As  initial 
conditions  for   (12.1)  we  have   therefore j 


(12.3) 


^l(O)   -i-ii   =^ 


w^(0)   s=  0 


and 


(12.4) 


W2(0)   ==  1  -  »3i 


^2(0)   =  w^^   , 
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In  which  w-^  and  Wg.  are  values  obtained  from  earlier  calcu- 
lations for  the  case  of  an  Initial  gas  pressure  ratio  of  125. 

The  accompanying  graph  shows  the  variation  of  w,  as  a 
function  of  the  time   r,  that  Is,  it  indicates  how  the  tube 
refills  after  the  valves  open,  i.e.  how  the  scavenging  phase 
takes  place.  The  numerical  integration  was  carried  out  up  to 
the  point  where  the  pressure  in  the  gas  layer  had  reached  its 
minimum*  The  value  of  w^  as  a  function  of  r  was  also  obtaiiiedf 
but  is  not  given  here. 

The  essential  conc!Luslon  to  be  drawn  from  the  graph  is 

that  the  refill  process  is  very  slow.  As  indicated  on  the 

graph,  the  time  of  stroke   t    required  for  the  gas  to  expand 

s 

from  a  pressure  p  *=  125  p^  down  to  p  =  p^  with  the  valves 

e  e 

closed  is   T  ST  0,205.  When  the  valves  are  opened  and  the 
refill  flow  (or  scavenging)  takes  place  we  observe  that  only 
1/5  of  the  exhaust  chamber  has  been  refilled  during  a  time 
of  about   7"=  0»6,  In  other  words,  the  scavenging  operation 
(as  we  have  considered  it  here)  would  require  three  times  as 
much  time  as  the  power  stroke  and  even  then  the  refilled 
volume  would  be  less  than  that  of  the  gas  layer  at  the  instant 
of  opening  the  valves* 

That  the  above  conclusion  was  to  be  expected  on  physical 
grounds  is  clear:   The  mass  of  water  in  front  of  the  valves 
is  at  rest  when  the  valves  are  open.   Once  the  valves  are 
opened,  this  mass  begins  to  move  under  the  pressure  difference 
existing  between  the  gas  pressure  and  that  in  the  water  in 
front  of  the  tube.  Initially  this  difference  is  very  small. 


A  few  calculations  were  made  in  other  cases  also,  but  we 
do  not  report  on  them  here  in  detail.   However,  one  result 
may  be  of  interest.   It  is  that  the  flow  into  the  exhaust 
chamber  is  not  greatly  affected  if  it  is  assumed  that  the 
valves  open  when  the  gas  pressure  has  dropped  to  a  value  of 
0.9  p^  Instead  of  p^. 
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so  that  the  velocity  of  the  mass  of  water  increases  ve}?y 
slowly  at  first  (see  the  attached  graph).   The  pressure 
difference  is  built  up  only  as  the  gas  layer  over-expands 
due  to  the  motion  of  the  water  in  the  exhaust  chamber,  and 
this  appears  to  be  a  rather  slow  process. 

It  seems  certain  that  efficiency  of  operation  of  the 
hydropulae  requires  that  the  scavenging  phase  should  not  take 
up  too  large  a  fraction  of  the  entire  cycle  of  operations. 
For  the  very  simple  design  of  a  hydropulae  discussed  above, 
it  is  then  clear  that  the  section  of  the  tube  in  front  of  the 
valves  should  be  eliminated  entirely  so  that  the  refill  would 
take  place  as  quickly  as  possible i  Clearly,  the  smaller  the 
initial  mass  in  front  of  the  valves  the  greater  the  acceleration 
at  the  beginning  of  the  process  and  the  more  rapid  the  entry 
of  water.  However,  if  the  front  section  of  the  tube  is  very 
short,  or  nonexistent,  our  assumption  that  the  pressure  over 
the  open  end  is  constant  is  likely  to  be  quite  wrong.  To 
obtain  reasonable  accuracy  in  such  cases  it  would  be  necessary 
to  consider  the  flow  pattern  near  the  end  of  the  tube  in  a 
detailed  way  —  both  theoretically  and  experimentally.  The 
conclusion  that  the  part  of  the  tube  in  front  of  the  valves 
should  be  eliminated  entirely  to  aid  in  efficient  scavenging 
is  just  the  reverse  of  the  conclusion  reached  in  the  preceding 
section,  where  the  operation  of  the  valve less  hydropulae  was 
considered.   There  we  found  that  a  long  constricted  front 
portion  of  the  tube  was  favorable  for  the  operation  of  the 
valveless  hydropulse,  principally  because  a  high  effective 
inertia  for  the  mass  of  water  in  the  front  part  of  the  tube 
would  be  obtained.  For  scavenging,  just  the  reverse  effect 
is  desired,  which  leads  us  to  suspect  that  the  front  section 
of  the  tube  ought  to  be  flared  out  to  promote  scavenging,  as 
Indicated  on  the  flgtiro: 
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-Valv 


Direction     of      motion 


Another  of  our  basic  assumptions  lihich  may  be   seriously 
at  fault  under  the   conditions  of  the   scavenging  process  is 
that  referring  to   the  existence  of  a  layer  of  gas  with  ylane 
interfaces  across   the  whole   section  of   the   tube.      Once   the  gas 
pressure   has  dropped  down   to  about   the  exit  pressure,    the 
actual  configuration  of  the  gas  volume  may  differ  radically 
from  that  of  a  slice.     It  seems  possible   that  the  gas  volume 
may  break  up  Into  bubbles  clinging  to   the   top  of   the   tube, 
and  that  gas  might   then  escape   through  the  front  end  of  the 
tube,  especially  if  the  valves  are  close   to  the  front.      This 
is  a  question  which  should  be  dealt  with  experimentally. 

Of  course,    the  very  rudimentary  design  upon  which  our 
above  remarks  are  based  is  not  the  only  possible  one.     In 
fact,   our  discussion  indicates   that  the  kind  of  "passive" 
scavenging  encountered  in  our   simple  model  should  be  replaced 
by  some   sort  of  "forced"   scavenging.      That  is,   part  of   the 
energy  in  the  gas  should  be  used  to  Increase   the  rate  of  flow 
during   scavenging.      The   turbo-jet  offers  an  example  of  a 
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solution  to   this  problem  for  duct  devices   in  air,   though 
here   the  purpose  of  the   turbine   is  more  to  create  a  high 
pressure   in   the  combustion  chamber  for  the   sake  of  thermo- 
dynamic efficiency  than  to  increase   the  flow  in  order  to 
increase   the   thrust -augmenting  mass-     Just  what  means  should 
be   taken  to  accomplish  similar  results  for  the  hydro-devices 
is  far  from  clear  at  the  present  stage  of  development  of   these 
devices.      Perhaps   valves  distributed  along    the   aides   of  the 
hydropulse   tube  might  be  used  to  facilitate    the  refill  process* 
Another  method  which  might  pi^ovide  forced  scavenging  is   the 
following*     Two   tubes  placed  aide  by  side  would  "fire"   into 
a  common  nozzle,   the   tubes   being  fired  alternately.      The 
high  velocity  and  consequent  low  pressure  produced  in  the 
nozzle  by  the  firing  of  one  tube  would   then  facilitate   the 
scavenging  of  the  other. 
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Section  15.  Effect  of  finely  distributed  gaa  bubbles. 

In  actual  practice  It  may  be  that  the  high  pressure 
gases  will  not  be  Introduced  as  a  single  layer  but  rather  In 
the  form  of  a  large  number  of  small  bubbles  thoroughly  mixed 
in  the  water.  For  example,  this  might  occur  if  th9  chemical 
which  reacts  with  the  water  to  form  the  gaa  were  to  be  thrown 
into  the  water  as  a  finely  ground  powder.  It  is  of  interest 
to  know  how  this  might  influence  the  behavior  of  the  hyd re- 
pulse, as  compared  with  the  behavior  when  the  gas  is  in  a 
single  layer. 

The  following  assumptions  will  be  made: 
1)  The  gas  bubbles  are  distributed  uniformly  in  the 
water. 


water. 


2)  The  gas  bubbles  move  together  with  the  surrounding 

3)  The  initial  velocity  of  the  water  in  the  tube  is 
zero. 

In  other  words,  the  mixture  of  gas  bubbles  and  water  will  be 
treated  as  homogeneous. 


Q 
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We  shall  ahow   that  suoh  a  mixture  of  gas  bubbles  and 
water  Is  equivalent  as  far  as   the  behavior  of  the  hydropulse 
is  concerned  to  a   single  gas   layer  with  a  different  total 
volume  of  the   tube*     More  precisely,   let  V™  be   the  voliame  of 
that  part  of  the  water  which  contains   the  gas  bubbles.      Then 
only  1/Z  of  V^  is  effective  as  far  as  thrust  augmentation  is 
ooneernedZ      That  is,   if  V_   Is   the  total  volume  of  water  in  the 
tube,    the  equivalent  tube  with  a  single  gas  layer  should  have 
the  volxjme 


Ve   -  §  V^  ' 


instead  of  V     -  V^.      The  reason  for   this  is   tlmt  the  mixture 
of  gas  bubbles  and  water  does  not  partake  of  the  full  velocity 
of  the   clear  water,   but  has   a  velocity  which  varies   linearly 
from  zero  at  the  front  end  to   the  velocity  of   the   interface. 
Of  course,    the  factor  1/3  is   to  be  considered  as  a  rough 
estimate,   since  any  deviations  from  the  assumptions  1),  2) 
and  3)  would   change  it. 

At  any  time   t,   let  Vq  be   the  total  volume  of  the  gas 
included  in  all  the  bubbles,   and  let  Vg  be   the  volxirae  of  the 
gas-water  mixture.     Since  the   total  amount  of  water  in  the 
mixture  remains  constant,  by  assumption  2),  w«    have   (see 
the  accompanying  figure): 

(13.1)  ^B  ^  ^0  *  V 

Furthermore,   the  density    Pg  of   the  mixture  is  very  closely 
given  by 

(13.2)  p^^  p  ^  , 
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in  lihlch    p   Is   the  density  of  the  liater* 

Let  us   compute   the  velocity  of  the  mixture  at  any  cross< 
section,  enclosing  with  the  front  end  the  volume  V.     At  the 
front  end,   the  velocity  of  the  mixture  is  0;     at  the   interface 
between  the  mixture  and  the  clear  water,    the  velocity  is 

• 
V 
0 

in  which  V^  is   the  time  rate  of  expansion  of  the   total  volume 
of  gas   (in  all  hubbies)  and  Q  is   the  cross-section  area  of  the 
tube,   assumed  constant.     Assumption  3)   is  used  at  this  point. 
We  proceed  to  show   that  the  velocity  of  the  mixture  varies 
linearly  with  V:      The  rate  of  increase  of  the  volume   of  the 
gas  enclosed  between  this  cross -section  and  the  front  end  is 

V    •  J. 

and  thus  is  equal  to  QU  where  U  is  the  velocity  of  the  mixture 
of  the  section  characterized  by  the  volume  V.   Thus 

(13.3)  TJ=^.^ 

which  varies  linearly  with  V. 

The  kinetic  energy  of  the  mixture  is 


f'  I  ^B^^v  .  \,\  .  fe)^    f  v^av  =  i  4  .  (§  V,). 


The  kinetic  energy  of  the  clear  water  is 
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The  total  kinetic  energy,  making  use  of  (13.1)  ia  therefore 

•  2 

(13.4)  IpJ-  (V^  -§V^  -  V(j). 

Expression  (13.4)  is  identical  with  the  kinetic  energy 

In  the  case  of  a  hydropulae  with  a  single  layer  of  gas  when 

the  initial  velocity  of  the  water  in  the  tube  is  zero,  except 

that  V  -  ^  V«,  replaces  V^.  Furthermore,  the  work  done  hy 
e   w  w  e 

the  gaa  is  the  same  in  the  two  cases  since  it  depends  only  on 
the  adiahatlc  law  and  the  change  in  volume  of  the  gaa.  The 
energy  equation  thus  shows  that  the  equivalent  hydropulse 
with  a  single  block  of  gas  has  a  volume 


^  2 
e 


v^  -  "S  v^- 
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Appendix  I 

The  method  of  finite  difference 3  applied  to 
the  nimerleal  solution  of  hydropulse  differential  equations* 

In  certain  parts  of  this  report,  (notably  in  Sections 
9,  10,  11,  and  12),  it  was  necessary  to  determine  numerically 
the  solutions  of  systems  of  non-linear  differential  equations 
under  prescidbed  initial  conditions.  Further  extensions  and 
refinements  of  the  hydropulse  theory  may  well  require  the 
integration  of  still  more  complicated  systems  than  those  we 
have  considered  hitherto.   Consequently,  it  i»as  thought  worth 
while  to  include  a  brief  discussion  of  the  methods  we  used 
and  found  convenient  for  this  type  of  problem,  even  though 
the  general  method  of  numerical  integration  by  using  finite 
differences  is  very  well  known. 

The  differential  equations  of  motion  which  appear  in 
Sections  9,  10,  and  12  have  the  general  form: 

^  -  f (u,v) 
(I.l)  dx*^ 

t^      XT 

w  -  g(u,v)  , 

dx'^ 
with  initial  conditions  (i.e.  conditions  at  x  =  x. ): 


(1.2) 


^U=x^  =  ^i      ^lx=x^  "  ^1 


duj     _  ,  »     dvi     _  „» 
3x'x=x,  ~  ^1     ^lx=x,  -  ^1 


^1    •"     ^  ^^1 
in  which  xx^,  m^,  v^,   and  v^  are  given  constants.  We  wish  to 
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determine  u  and  v  as  functions  of  x*^  It  is  possible  to  solve 
equations  (I.l)  numerically  by  replacing  them  with  appropriate 
difference  equations.   This  can  be  carried  out  in  the  following 
way; 

Let  us  consider  one  of  the  dependent  variables,  u« 
The  value  of  u  at  x  +  h  and  x  -  h  can  be  expressed  in  terms 
of  the  value  of  u  and  its  derivatives  at  x  =  x^  by  a  Taylor 
series. 

(1. 3a)     u(x^j+  h)  =u(x^)  +hu'(x^j)  +  |h^u  (x^^) 

1  ^3  «  ,_  >  .  1  v,4  IV, 


+  ^  h-^u  (x^)  +  ^   h'u-^Mx^)  +  ... 


(1. 3b)     u(Xq-  h)  =  u(Xjj)  -  hu*(x^)  +  ^h\   (x^) 

Adding  (I. 3a)  and  (I. 3b)  we  have 

u(x^j+  h)  +  u(Xq-  h)  -  2uCx^)  +  liV(Xjj)  +  ^^^^^ ^^q)   +  ••• 

or 

o  w       1   4  TV 

(1.4)  u(x^^+  h)  =  2u(x^)-  u(x^-  h)  +  h'^u  (x^)+^  u  (x^)+  ... 

If  h  is   taken  small  enough.    It  is  possible   to  neglect  all 
terms  of  order  higher  than  the   second  In  h,  yielding 

(1.5)  u(x^+  h)  =  2u(x^)   -  u(x^j-  h)   +  nVd^) 

If  the  values  of  u  and  v  are  both  known  for  i  *=  x^,   than  W6 
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n 
can  use  equation  (1,1)  to  obtain  u  (x  ): 

u"u^)  =  f(u[x^],  v[x^])  =  f(x^) 
and  hence  (I«5)  can  be  written  In  the  form 

(1.6)  u(x^+  h)   =  2u(x   )   -  u(x  -  h)   +  h^f(z    )   , 


with  t(x^)  known. 

If  we  take  x^  =  x^  in  (I. 3a)  we  can  determine  u{x.+  h) 
and,  similarly,  v{x^+  h)  with  the  help  of  equations  (1,1)  and 
(1.2)  thus 


u(x^+  h)   =  u^+  ^1  +  ^  h^f(u^,v^)   , 
(1.7)     and 

r{x^+  h)  ^  y^+  hy[  +  ^  h^g{u^,^r^)   , 

in  which  derivatives  of  03?der  higher  than  the  second  are 
ignored.  It  is  now  possible  to  determine  u(x^+  2h)  and 
v(x^+  2h)  from  equation  (1.6)  and  the  corresponding  equation 
for  V,  by  taking  x^  =  x^  +  h.  Taking  x^  in  succession  equal 
to  x^  +  2h,  x^  +  3h,...  ,  x^  +  nh,...  it  is  possible  in  this 
way  to  calculate  approximate  values  for  u  and  v  as  functions 
of  X  at  points  which  are  a  distance  h  apart. 

Any  integral  of  a  function  of  u  or  v,   say 

,b 
I  =     I    d)(u)dx 


=     1    <^(u)d: 


can  be  determined  quite   simply  with  the  help  of  Simpson's 

One  could  include  higher  order  tenns.     For  this  it  would 
only  be  necessary  to  calculate  higher  derivatives  of 
u  at  X  =  Xjj  by  differentiating   (I.l). 
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rule,  ualng  intervals  vhich  are  some  integral  multiple  of 

h  in  length,  since  the  values  of  u,  and  hence  of  <p  {n),  are 

known  at  the  ends  of  all  such  intervals. 

In  order  to  expedite  the  integration  of  (I.l)  it  is 

desirable  to  use  the  largest  value  of  h  consistent  with  the 

accuracy  which  is  demanded  in  any  particular  instance*  The 

largest  term  in  equation  (I»4)  which  is  neglected  is  ygh  u  (x  ) 

and  this  will  be  the  approximate  error  in  any  one  step» 

This  error  can  be  estimated  in  a  simple  manner  in  the  course 

of  the  numerical  calculations  from  the  following  considerations. 

n 

Equations  (I.3),  (1.4)  and  (1.5)  can  be  written  for  u  rather 
than  u  as  the  function  of  x.  We  get 

(1.5*)   u"(x„+  h)  =  2u''(x„)  -  u"(x„-  h)  +  hV^(x^) 


or 

(1.8)  ^hV^(x^)  -  ^  li^u"(x^+  h)  -  2hV(x^)  +hV(x^-  h) 

=  ^  h2f(x^+  h)-2h^f(xQ)+h^f(x^,-  h)  , 

which  is  ^  of  the  fourth  difference  of  u.   Thus,  by  taking  a 
running  second  difference  of  h  f(x)  we  obtain  an  estimate  of 
the  error  in  any  step.  If  this  error  is  much  smaller  than 
necessary  we  may  increase  h  and  carry  out  the  integration  more 
rapidly.  If  the  error  is  too  large,  h  should  be  decreased. 
Any  such  change  in  h  from  h^  to  hp  during  the  calcu- 
lations, say  at  X  =  a,  involves  an  interpolation  procedure 
to  determine  u  at  x  =  a  -  hg  before  (1.6)  can  be  used  again. 
This  interpolation  can  be  carried  out  in  several  ways,  for 
example,  by  the  use  of  (I. 5b)  with  suitable  estimates  for 
u  (a),  etc.  (it  should  be  noted  that  the  integration  of  u 
and  V  must  be  carried  out  simultaneously,  and  the  same  h  should 
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be  used  for  both  Integrations*) 

The  follovlng  numerleal  example  vlll  Illustrate  the 
nay  in  which  the  choice  of  h  affects  results! 

We  consider  the  case  of  tvo  gas  layers ,  as  described 
in  Section  10.  The  equations  of  motion  for  the  system  were 
given  as 


(1.9) 


^8^"^  -fs'^"^  -2a^" 


rl^J)''^   -1=  2{<f"  +77")(1  -<f  -Voc) 


or 


(I.IO) 


^"=^(^1^"''-^^^) 


a 


''"=^rT=?W'^.''-^-^'-^'^>-''-''>-"'^ 


We  consider  the  case  in  which  the  parameters  have  the 
following  values; 


(I.ll) 


7=  1.5  ,   ^^  «  77g  -a  =  0,25  , 


8 


and  as  initial  conditions,  «e   take 


(1.12)        ^^  =  0.01,     r)^  5=  0.025,      S,^  -     n[  =  0 


for      7-=  0 


Substituting  in  (1.9)  we  obtain 


(1.13)  ^^  =  186.76,     Tl^  -  "  165.35 
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and    (I.V)   becomes 
(1.14) 


^(Ar)   =  .01   +  93.58(At)2  _  3.1795  x  IO^(At)*  +  .,, 


77  (At)   =  ,025  -  82.68(Ar)2   +  3.2712  x  lO^(Ar)*   +  ... 

In  order  to  make   these   series   converge  fairly  quickly, 

we  must  take  At    of   the  order  of  0.005  or  less.      Once  At    has 

heen  chosen,   we   can  proceed  as  already  indicated,  with  the 

repeated  use   of  equation   (1.6)   in  which  u  and  v  are   to  be 

replaced  by    <^    and      77  . 

The  accompanying   table   shows  the  results  which  are 

obtained  for    <f    as  a  function  of   t    when  h  (or  At)   takes  on 

the   values   0.00125,   0.00250,   0.0050,   0.0100. 

It  can  be  seen   that  h  =  0.00125  is   smaller  than  necessary 

for  our  purposes,   and   the  accuracy  is  not  appreciably  affected 

by  taking,   rather,  h  =  0.0025.     It  will  be  found,  furthermore, 

that  ^  of   the   fourth  difference    is  completely  negligible   in 

the  former  case    (  At    =  0.00125)  and  even  over  most  of  the 

range  for  h  =  0.0025. 

In  order  to  illustrate    the  effect  of   changing 

the   interval,  we  have   carried   this  out  for   the  case   in 

which  h  =  0.0025  initially.      The   interval  is  successively 

doubled   at      7"  s=  0.0350  and  at      Tri  0.050,   as   the  fourth 

difference   is  quite   small  at  each  of   these   points.      It 

can  be  seen   that  the  error  introduced  as   a  result  is  quite 

small  —  of   the  order  of  1  per  cent.     It  should  be  noted, 

in  connection  with  the  use  of  equation   (1.6)  at 

T    =  0.0350,   that  when  the  AT  Interval  is  doubled  we 
o  ' 

must  take 

^(  :r    -  h)  =  0.02870 
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^ 

^ 

- 

r 

Ar  =0.00125 

Ar=0, 

►  0025 

At  =0.005 

At  =0.01 

0. 

.01000 

0.00125 

1015 

.00250 

1058 

.01057 

.00375 

1127 

.00500 

1218 

1217 

1214 

.00625 

1327 

•  00750 

1450 

1449 

.00875 

1583 

.01000 

1722 

1720 

1718 

.01616 

.0125 

2001 

.0150 

2266 

2261 

.0175 

2494 

.0200 

2670 

2657 

2565 

.0225 

2788 

.0250 

2853 

2827 

2718 

.0275 

2876 

.0300 

2870 

2837 

.0325 

2846 

.0350 
.0375 

2814 

2782 

2781 

.0400 

2753 

02756 

2729 

2674 

.0425 

2735 

.0450 

2731 

2738 

2719 

.0475 

2743 

.0500 

2773 

2784 

2775 

2809 

.0525 

- 

2822 

.0550 

2890 

2907 

.0575 

2977 

.0600 

3083 

3113 
(Ar=0.01) 

3115 

3222 
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eorres pond leg  to  h  ~  0*005,   Instead  of 

^(t^  -  h)  -  0.02846 

uhlch  would  have  been  correct  if  h  had  remained  unchanged. 

From  Inspection  of  the  table.   It  seems  reasonable   to 
conclude   that  It  would  be  most  efficient  to  choose  h  «=  0.005 
from  the  outset  If  an  accuracy  of  about  1  per  cent  Is  desired. 
Similar  considerations  would  determine   the  most  efficient 
choice  of  h  under  different  Initial  conditions. 
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Appendix  II 
Special  method  for  oalculatloa  of  y(t)« 

The  method  of  finite  differences  outlined  in  Appendix  I 
above  can  of  ooiirse  he  used  to   treat  simpler  cases   than  those 
of  Sections  9^   10,   and  12*     In  fact,  «e   have  applied  it  to  find 
solutions   in  the   simple  case   of  the   straight  tube  with  one   layer 
of  gas*     However,  for  these   cases  the  difference  method  is  more 
laborious   and  time-consuming  than  necessary*     In  this  Appendix 
we  outline  a  method  of  calculation  for  the  case  of   the  straight 
tube  with  one   layer  of  gas  which  makes   it  possible   to  determine 
accurately  the  gas  volume  as  a  function  of  the  time  over  the 
whole  stroke   in  a  fraction  of  the   time  required  by  the  method  of 
finite  differences*      That  the  method  is  very  accurate   (with  an 
error  less   than  1  per  cent  in  general)  we  know  by  comparison 
with  the  solutions  by  finite  differences* 

The  differential  equation  to  be   integrated  la  equation 
(4*9)  of  Section  4  for  the  case  q  =  1,   as  follows i 

(II«1)  3^(v  ) 1-^^  , 

t   dv 

in  which  V  =  ^«  the  variables  being  in  the  dimensionless  form 

discussed  in  Section  4*  Since  g=(v  )  =  2v  we  may  write  (11*1) 
in  the  form 

«   v>-">'-l 

The  differential  equation   (11*2)  is  to  be  integrated  under  the 
initial  conditions t 


(II.3) 


v(0)   «=  Vj^  , 
T  (0)  F  r^ 
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By  integration  of  (II. 1)  we  obtain  v  aa  a  function 


of  v: 

(II.4)         v*2  ^  v^^  + 


t 


i 


in  which  Vm   is   the   initial  value  of  ▼  ^  as  ^oove. 
We  note  that 

It  follows   that  all  derivatives  of      t(v)  with  respect  to  v 

can  be  obtained  by  making  use  of   (11*2)   since  v     can  be 

considered  a  known  function   of  v  by   (II»4).     In  what  follows 

it  is  convenient  for  various  reasons   to  determine    t    as  a  function 

t         « 
of  Y  rather  than  v  as  a  function  of      r.     By      "^   »     '^   »  etc. 

we  mean  differentiations  with  respect  to  v« 

The  method  of  approximating   the  function    t  (v)  which 

we  propose  is  based  on  our  knowledge  of  the  general  character 

of  the  function:     We  know  that    r(v)   is  a  nearly  linear  function 

of  V  over  an  interval  v,  +  e  <▼<▼-#  ^>  0  which  extends  up  to 

v_    (the  gas  volume   corresponding  to  the  exit  pressure  p   )   but 

not  all  the  way  down  to  the   initial  value  v-   of  v.     Near  v^^, 

in  fact,   the  c\jrve    t(v)  has  high  curvature.     Our  method  of 

approximation  thus  consists  of  using  a  series  development  for 

T  In  powers  of  v  -  ▼     for  all  values   of  v  in  an   interval 

v^  <  v-  <  V  <  T  ,   in  which  v-  is   to  be  suitably  determined, 

together  with  a  different  type  of  approximation  valid  over  the 

remaining  interval  from  y.    to  v,»     The  main  difficulty  arises 

in  piecing   these  two  different  types  of  representation  of  the 
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solution  together  at  v  -  v,.  The  circumstances  are  Indicated 
in  the  acconpanying  figure,  which  should  he  compared  with  the 
accurate  solutions  given  by  Figure  (7»10)  (cf«  the  later  parts 
of  Section  7)j 


To  carry  out  this  program  we  begin  by  expanding  t(v) 


about  ▼  =  V  : 


(II.6)   r=  T^+  ^i(-V3)  +^(v-Y3)2+^(v.Vg)3+  ,.., 


with 


(II.7) 


a 


r  9 


lA, 


i^.U".'  *  wT? 
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Taj  use  of  (II. 5)  and  (II»2),  The  quantity  r^  ia,  of  courae, 
th©  "time"  of  stroke.  The  quantity  ▼.  is  Imo^n,  but  t  and 

t  I        B  8 

T     are  not  known*      However,     r     can  be   computed  using   (II.4), 
since  V     =  1/^  ~  1/^5     '®   need  only  expand   (1-^  )"     in  powers 

of    ^    and   then  cbttj  out  the  integration  in   (II*4}  to  obtain 

to 
V       with  any  desired  accuracy  as   a  power  series   in  v.      Once 

this  is  done    the  coefficients  in   the  series    (II. 6)  for    r   would 

all  be  known  with  the  exception  of  the  first  one     t  .      To  fix 

the  value  of      r    we  need   to  know  only  one  pair  of  values    (v,  r), 

B 

with  V  a  value  not  too  near  v^  (or,  r  not  too  near  zero)  so 
that  (II. 6)  still  converges  reasonably  well. 

We  next  investigate  the  behavior  of  t(v)  near  v.. 
Prom  T  =  -j  7  n  and  (II.4)  we  see  that  the  behavior  of  r(v) 

»     t 
near  v  ~  v,  depends  strongly  on  the  initial  value  v.   of  v  8 

I     ^       t  1  •"■ 

If  v^  ~  0,  then  r  behaves  like  — —  when  v  is  near  v. ,  but 

if  V,  ?^  0  then  r  is  regular  for  v  =  v..   Thus  ''"  would  possess 
two  different  developments  near  v  =  v*,  as  follows: 

i 

(II. 8)  For  v^=  0:  T=  a  (v-v^)  +b(v-Vj^)+  c(v-Vj^)  /^  +  .«• 


(II.9)  For  v^?^  0:  t= — ^  -f  j3{vY^)^+    (f(v-Vj^)^+  e  (v-v^)*+  ... 


If  we  were  to  compute  the  coefficients  a,b,c,/S,d  ,  ^  by 
evaluating  the  successive  derivatives  of  '''  at  v  =  ▼,  ,  we  would 
find  the  series  to  converge  very  badly.  Consequently,  we  break 
off  the  developments  with  the  terms  in  (v-Vj^)  /  and  (v-v^) 
and  compute  the  coefficients  in  the  following  way:  The 
coefficients  a  and  y8  are  determined  so  that  t  has  the  correct 
value  at  ▼  ;=  v^,  but  b  and  c,  d  and  €   are  chosen  in  such  a  way 
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that  the  first  derlvatlya       r     has   eorz>ect  yalnes  at  each  of 
the   two  points  ▼  -  ▼«  ^^^  "^  '^  ^3  gl"^®°  "^ 

(II.IO)  ^2  "  ^   -^  ^^  •   ^3  =  ^1  -^  ^^  • 

In  this  way  the  point  v-  referred  to  above  is  located*   The 
values  of  r    (vg)  and   r  (v,)  are  calculated  from  (II«4)  and 
the  relation   t  s  1^   , 

The  coefficients  a,b,c  can  now  be  calculated  as  follows  J 
One  finds  readily  that 

"    2 
Lim  V  =  -w  , 

a** 
v-^  V- 


so   that  a  =J2/Vt    •     Hence  a  is  fixed   since  v.   is  known. 


Inserting  this  value  for  a  we   then  obtain   the  following  pair 
of  linear  equations  for  the  determination  of  b  and  cs 


^•{vg)  =  ^JK  (y^-v,)-  2  +  b  *  I  c(v2-T,)5 


t'IVj)    =  I  /X  (T3-V^)-  I   +  b    +  I    C(T3-T^)?    . 


The  calculation  of    ft   ,    ti    ,  and    e    follows  the  same 
scheme* 

The  formulas   (II.8)  and   (II»9)  are  found   to  be  quite 
accurate   approximations   to     r  (v)   in   the   interval  v.    <  v  <  v_» 
By  substituting   (Vj,  't^)  in    (II. 6)  an  accurate  value  for       r^ 


It  is   clear  that  Vg  and  v,  need  not  be   chosen  exactly  in 

this  way.     We  found  that  these   choices  for  Vg  and  v,  gave 
good  results. 
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can  be  found  —  aince  {II»6)  converges  quite  well  for  v  i"  v_« 

o 

Once   the  value  of        t    ±3  known,    (II. 6)  In   the   interval 

3 

V,  <  v  <  V  together  with  (II.8)  or  (II.9)  in  the  interval 

V.  <  V  <  V-  constitute  very  good  approximations  to   r  (v)  over 

the  whole  range  v,  <  v  <  v_. 

The  method  has  been  checked  in  a  considerable  number  of 

cases  by  comparison  with  the  results  obtained  by  finite 

differences,  with  results  that  indicate  the  error  in   t(v)  to 

be  less  than  1  per  cent  of   r  .   The  method  of  approximation 

s 

outlined  here    can   be   carried   through  by  an   experienced 
calculator  in   a   specific   case   in   less   than   one-half  hour,   while 
the  difference  method   requires   about  four  hours. 
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